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PREFACE 


An attempt has been made to keep the programs as subroutine oriented 
as possible. Usually only the main programs are directly concerned with 
the problem of total cross sections. In particular the subroutines 
POLFIT (IV. A), BILINR (V.B), GASS59/MAXLIK (VI. D) , SYMOR (VI. E), MATIN 
(VI. F), STUDNT (VI. G), DNTERP (VII. E) , DIFTAB (VII. F), FORD IF (VII oG), 
EPSALG (VII. H), REGFAL (VII. I), AND ADSIMP (VII.J), are completely 
general, and are concerned only with the problems of numerical analysis 
and statistics. Each subroutine is independently documented. 



i 


TABLE OF CONTENTS 

PAGE 

I. Discussion 

II. Uhits Il-l 

III. XSECT - Generation of Total Cross Sections . . m-l 

IV. PHASHT - Fitting of Glory Extrema . . IV-1 

A. POLFIT - Forward Curvilinear Regression • . IV-3 

B. Examples IV-5 

V. XTREMA - Fitting of Glory Extrema ..... V-l 

A. Fit V-2 

B. BILINR - Bilinear Regression. V-3 

C. EXAMPLES V-4 

VI. QFIT - Fitting of the Total Gross Sections VI-1 

A. TRIAL VI-4 

B. FINAL VI-5 

G. XSECT/DXSECT ; VI-6 

D. GASS59/MAXLIK - Nonlinear Regression VI-7 

E. SYMQR - Eigenvalues and Eigenvectors of a Symmetric 

Matrix. VI-9 

F. MATIN - Gaussian Elimination VI-10 

G. STUDNT (Student - t Statistic) VI-11 

H. Examples. . ; VI-12 

VII. JWKB - Cross Section Program VII-1 

A. PHASHT - JWKB Phase Shifts. VII-6 

B. VEFF/F/F2/F3. * VII-8 

C. POTIN/POT - Potential Routine ....... VII-9 



ii 

PAGE 

oo 

D. QJBVAL - Jeffreys-Born cLJr' ...... VII-12 

E. DNTERP - Finite-Difference Interpolation VII-13 

F. DIFTAB - Difference Table Construction.. VII-14 

G. FORDIF - Forward Difference Calculation VII-15 

H. EPSALG - e-algorithm Extrapolation VII-16 

I. REGFAL - Root of Transcendental Equation VII-18 

J. ADSIMP - Adaptive Simpson’s Rule Integration VII-20 

K. Sample Output ... VII-21 



1-1 


I. DISCUSSION 


The total cross section Q(v) , where v is the relative velocity, 
is assumed to be a sum of two components: 5 dependent upon the 

I 

long-range force constant C,. , and AQ (v) , an oscillatory term due 

o g 

to the presence of glory extrema 

Qtir) - Q S uL Clr) +• AQ^cv) 


(i-D 


where 


/ Ca \ % ls 

Q sll c^) ^ 8,083 ( — ) 

AQ, tv) - ±^ lz b ^ ( 27 „- fj 


(I-2a) 


k 2 - w- 2 '’) 

where E m JiL'V' /"^ is the wavenumber, L g is value of the orbital 


quantum number corresponding to > a maximum in the phase 

U 

shift, and is the second-derivative of with respect 

to evaluated at L . 

It can be shown that by expanding about the high-velocity limit, 
(2b) may be written 

AGM-v) = (w^) 5,l 3„ (^) ,2 [i+ H ^+ l 4 £ +"'] /uil ^ I - ; 


-3) 


r— — — 2. 

where E “ 2 . ^ r " is the energy, 


4 = 


3r +46^ (a< + 

4 -ft or 


Ai & ^ A z £ 
E E 2 - 




(1-4) 


For a more complete description of the derivations of the formulae see 
Wise. Theor. Chem. Institute Report WIS-TCI-469 (1972). 
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and r^ is the position of the minimum of (of depth e) the potential. 

The constants g Q , a^ are related to the curvature of the potential well 

at r = r , and the constants H. , H„ , ... , A. , A ', ... , 

m 12 12 

were found to be to a good approximation independent of the curvature, 

e , and r^ , and of the form of the potential. Furthermore, for the 

* * 

energy range 1 <_ E _< 100 , where E = E/ £ is the reduced energy, it 
was found only the terms and A^ , A^ , A^ and A^ were required. 

For a Lennard-Jones 12-6 potential, i. e . , 

Vtr) - e [ ?)‘] 

the values of the constants are: 


= 

0.186299 

(I-5a) 

o 



H l = 

3.267 +0.80 

(I-5b) 

a l = 

0.4215587 

(I-5c) 

u 

1 — 1 
<3 

- 0.1655 + .0016 

(I-5d) 

A 2 " 

0.1057 + 0.0160 

(I-5e) 

A 3 = 

- 0.0544 + 0.0160 

(I-5f) 

> 

•c- 

ii 

0.0139 + 0.0084 

(I-5g) 


In terms of free parameters, (1) may be rewritten 

Q Cor) zr Y, v~ ^ + Y L u /z (\+ ^ )^> <$> 

B 1 d-6) 

and 

A> , - UL + *j_ + Ik. ( A + Aze jJhll'kl*) 

4 or vE V 1 £ ^3 / (i-7) 
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/ C & \^ 5 

Y, = S,o83 ( -£) 

(2 W r m ?k Jo 

•» 

X 3 - 4- a ,er m /-ft 

Yu - 4 G*Y» 

^ (1-8) 


It has been shown previously that when the r term in the potential 
begins to affect the cross section, it may be accounted for by replacing 

Q sll ( v ) in C 1 ) b y Q SUj ( v ) + 2Tr/,S ® » where b “ c g / c 6 • 

It is clear from the form of (1-6) that Q(v) has an oscillatory 
component. The condition for extrema in Q is that sin (J) = + 1 , 
which can be shown to lead to: 

= z ^_( a 4 Qit + hi* 1 ' e \ A h^J 

&' If ^ t t g £7. J (1-9) 



where the maxima and minima in Q(v) are indexed by N(v ) in the order 

g 

they occur, i. e . , N = 1 corresponds to the first maximum, N = 1.5 cor- 


l 

responds to the first minimum, N = 2 is the second maximum, etc. The 

quantities are the associated values of v at the extrema. 
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II. UNITS 


All of the results and .programs below assume the following system 
of units: 

-24 

Mass M = 10 g = 1 ppg (picopicogram) 

Length L = 10 ^ cm = 1 A (Sngstrom) 

-13 

Time T = 10 sec = 1 dps (decipicosecond) 

In this sytem, 

1 amu = 1.660531 ppg 

1 a. u. = 0.52917715 $ 
o 

1 Km/sec = 1 A/dps 
1 eV = 160.21917 cpe 
1 Kcal/mole = 6.94793 cpe 

h = 6.626196 cpe-dps 
= 1.0545919 cpe-dps 
kT = 1.380622 cpe @ T = 100°K 
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III. XSECT 


The purpose of the main program XSECT is to calculate total cross 
sections Q(v) by means of (I — 6)-(I — 8). The assumption is made 
that and , A^ are those for the L-J (12-6) potential 

given in (1—5) . 

The first card of the data set contains (18A4). 72 columns of title 
information, which is reproduced at the top of the output. 

The second card of the data specifies N, MU, EPS, Al, BO, Cl, C6, 
in that order, according to the format (15, 6F10.4, E15.6), where: 

N - No. of velocities for which Q(v) is to be calculated 
MU - y (reduced mass) 


EPS - £ 


RM - r 

m 

Al - a (I-5c) 

BO, Cl - b Q and c^ are related to the leading terms in the expansion 


in powers of 1/E of L and 


in " 

8 ( m 


respectively. 
1/2 


Since the only quantity required is g Q = b Q /c^ , b Q 

may be set equal to g Q and c^ to 1 . 

C6 = C the long-range force constant, 
o 


The succeeding cards of the data contain the velocities V(I), 1=1 

N, for which the Q(v) are to be calculated, according to the format 
(8F10.4). 

The output of the program is a table of the values of V(I), QSLL, 
DELQ, Q, and QF for I = 1, ...» N, where 
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QSLL - Q SIi ( v) 

DELQ = AQ (v) 

§ 


Q ' Q(v > ■ Q S ll ( v > + ^Q g W 

QF ' A V v>/Q su. (v > 



XSECT 7 ; • •" - - • •- 

0 COMPILED BY FOR-V S7E6 ON 06 SEP 71 AT 17:4i:il. 

GE USED: CODE(l) 000270; D AT A ( 0 ) 000515 


PROGRAM TO GENERATE RAPID CROSS SECTIONS m 

R E AL ~ M U ~ r ' : • : pi ~~ ~ 

DIMENSION VI 200) » A (4) *TITL EI1B I M 

DATA A/-0.1S55* 0 .1 0 57 # -0 . 05 4 4 » 0 *01 3 9/ 1 Hl/3 .2 67/|7 . 

DATA PI/3 « 141 5926 5/ »HCR OSS /ji •054 5919/ 

READ 15*5) ITITLEI II *1=1*18) ' * 

FORMAT 1 1 8 A4 ) . 

~~ WR IT El'S *S1' _ TTITL E'HT : * X-r*T8T ' “ 

FORMAT ClHl»10X*ieA4l k- 

READ 15*10) N»MU»EPS»RM»A1*B0»C1*C6 

‘ FORMAT "ll5f6Rl0.4»ElS. 6 ) ’ $fS 

WRITE 16*20) N*MU»EPS*RM» AltBO* Cl * C6 

FORMAT I1H0*10X*»N0. OF VELOCITIES ** I8/11X *• REDUCED MASS** 
t~tl X *;* WELL-OEPT H-EPS** El 5 .6/1 IX* * POTTNTIAC 7TXNIMqTI"RM *TET5 

2 ’PHASE SHIFT COEFFICIENT A1 * * El 5.6/ 1 IX t ’BETA COEFFICIENT 

3 E15.6/11X* * COEFFICIENT Cl * » E15 . 6/1 1 X* * LONG RANGE C6**E15. 

“ READ I 5 *30 I T V II nT=I *R7 ~~:y. ’ “ 

FORMAT I 8 F 1 0 • 4 ) 

01 = C6/HCR0SS 

~D2 -■=—■ 2 « *Pr*RM*BO~ : |pn 

D3 r 4.*RM/HCR0SS |p 

D4 = 2.*PI*RM*HCR0SS/(EPS*C1I 

"WRITE (6*40)* T ~ 

FORMAT I // 1H0 » 5 X t * NO » * * 8X* • V • *1 3X * * OSLL • • 10 X * *D|L0 •• 11X • • 

1 10X**DELQ/GSLL*/4X.100C1H*I/1 J ’ 

D 0 2 0 0 I "= _ 1 * N T 

QSLL = 8o083*(Dl/V(I) l**C.4 / 

ETR r 0.5*MU*V(I)*V( I) 

~Z -eps/etr ~ ~ ' 

PHI = - 0* 75 *PI ♦ D3*EPS»IA1 ♦ Z*(A(1) ♦ Z*<A<2) ♦ ZMAI3I 

i Ai 4 m y n/vi i) 

“DELT1 -0Z*3QRT(04*Vtm* Cl* ♦ HI* Zl *3INI PHI1 

Q = QSLL ♦ DELO ’ 

QF = DEL Q/QSLL 

write f 6. SOT”I rV'l iTtm L'lD ELffWSr” 

FORMAT I 3X*I5*5E15 *6 ) 

GO TO 100 

-END 1 


El 5.6 


BO* * 


QTOT’i 


♦ Z* 


ND OF compilation: 


NO DIAGNOSTICS, 


compilation time; 

I7 0REQ OESTS: 

i/o words transfered: 


0.653 SEC, 
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IV. PHASHT 


The purpose of this program is to produce a fit of the experimental 
3 

variable Y = (N - -jr)v„ to the model 

o N 


Y = i + 


+ ••• 


e e 2 - 

where, after regression, the following ansatz is made: 


(IV-1) 


y _ 2a, 

ith 


j S- L - ' L ± ^ A; , 

O J / 


/ TT 'ft 


(IV-2) 


which leads to 


a, € ^ 


_ 7T h -£ 

2 . 

e' 2 -^ = 'TTfc S,/A, 

2. 

A* Si 


e - 


^2. Sj 


(IV-3) 


Values of £ r and £ are calculated from (IV-2) assuming A,, A„ , 
m 12 

A^» A^ are given by (1-5). "ERROR" denotes 95% confidence interval 
halfwidths. The confidence intervals were calculated using first-order 
propagation of error, and took the error in A^, A^ , A^> A^ into account. 


** 


I. e. , if 


^ - -f Cx, , then 


- (E,) 2 < + 1 + 
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Input to the program is as follows: 

i 

Card 1: title information in columns 1-72 (18A4) 

Card 2: NOB, IWGHT , MU (215, F10.4) 

Card 3 (or more): V(I) , I = 1, ... , NOB (8F10.4) 

' Card 4 (or more): INDEX(I) ,' I = 1, ..., NOB (8F10.4) 
where 

NOB - number of observed extrema 

IWGHT - if 0 data has assumed constant absolute error. 

data is assumed to have constant relative error. 
MU = y (reduced mass) 

V(I) - extrema velocity v^ for I-th extremum 
INDEX(I) - extrema index N for I-th extremum 


If 1 , 



R*SI PMASHT 1 ' 

tfc ooo toHriufeb #r fo«-v 


ON Or JAN It At 061*31*** 



DIMENSION XI 30) iT(30)*R 

<3O)»Vi30TtAlR>*^ 

DATA A ) /0* 12^.550/ *A/-0* 

U5s*Oti0S7i*biO9A 

/ ,EWr/0. TSFsit -Ay ,£ R a 1 A 

2/151031*27 

HE AD | N EXPERJ n£NT T j TC 


READ (5«5) 1 TITLE i|) #*• 

*«»■» - s 



inn 


WRITE (A t *l lT|TUUl»|*tl|6t V 

F ORMAt 1 1 H ! t ) OX » t OA *t> )., ;.V.v 

read- jn(no. of da n m nttvort w*#t~ , ... , . 

e rror Ino wi i > » v -;4f *< ti • t/*l IMP* n 
And reduced mass nd« •“* „.r '- 


NOTE #•* IF V IN UNTT-S Of *M/StC AND HU |N fl .C-fN ORAM, 

OF ENEROt «iU BE i*E-|R Er 6$* (NOTE THAT 0ON5 HAS UNITS Of ACTIO 
I*£-27 ERO-SEcI ■ t 

^XAV T$iiV 7 NO eTnJfiMTTMU ~~ - ' ■ ■ J~" “* .. , v; -y-vT^y 

FOR HAT ( 2lS*F IO*R I : *'■ ; •" v 

WRITE (6>20 ) NqB» (RGHT jMu 


FONWAT TiHQ*IOX**NO* OF OBSERVATIONS **•» 

ih *•••»! to/i i x i + reduced mass ...*»c)s*4)„ 

READ IN EXTREMA VELOCITIES AND ExTrEHA |N»'lCE$ 


li ’ iH'in.m.iAniiiCLsyii 



^I^Si^FORMAt, ( ////SxY* NO* * » «X . * V * » l §* ♦ ♦ ETR»* )| 
■ V|/AX#ao( IMF) /I • V: 




ft*® 


T U I * t INDEX ( 1 1 - 0*375) *V 1 1 ) 

'SCI I . • :!* /-. • • 

v I f MR«wiT *eo. *1 «< I ) • Mtri I r*y 1 1 ) ) 


i 

s 

« 


WRITE (4*401 NpEG *S I 0 ... -v . • 

1 . r0«HA:|^i/'/i0 ; x • ♦ • 1 DETCRHiNEDVioxnRMS * 

iSiOUAL ; I V« * . fct ***** I Its * 1 ***UlI**H : 


L LUCIi 


ENT » i i 0* I • ERROR MSI » S6 


SUITE; (0*70) J.C(||»EU> 






if ' : : /v ^v; i :*.•-• . , .. >.•■ ■. ; -- 

CM I fglrf ;• ■•••• ■ >T : . - -"x 


FORMAT i?X, UtElSiS »5X V C 1 5*5) ■ 

A£R i CdNSimi” : v? ~ 

EE * CONS*E ( 1 1 • 

EFSRH • AER/A 


EE2« EE/A 1 
WRITE (5*601 AERIES* EPSON *EE2 

FORHAt t/10X»*Al*EPS«RM ■ * , E | 2$ ,2 X . M ERR OR* • , E 1 2 . 5 . » )*/)0X*«CP! 


OM Al • Q«R2|65V .*.*,El2 t B*2X*MERR0f •« «C 1 2* 5 * * ) • ) 








> *4 *4 

>4 6 

5 3 • 
54* 


IF (NDEG • EQ • 0> GO TO lQO 
EER a CONS*C l 2 ) / A ( 1 ) 

— L 

I 

>47 

>50 

>51 

55* 

56* 

■ 57 \* 


EE » ABSfC0NS*E<2)/A( 1 I I 

EPS B EER/EPSRM 

EE 2 8 £PS*SQRT(RU . 1 1/C* l >**2 

♦ R(2.21/C12»|*2 -2 » *R ( I , 2 I /( C ( I ) *c 

>51 

58*f 

- 

12 U ♦ E«Al) 

•*« j 

>52 

59* 


W RITE { 6 s 85 ) E£R,EE»EPS»EE2 

it • • ■ 

>60 

60* 

85 

FORMAT ( /10X . ’EPS* EPS *RM = * ,£12 

• 5.2X » * (ERROR! *«.E12.S,» >*/|OX,«Ep 

>60 

61 * 


1 FROM Ai *Q»42lS59 '.«••• E.l 2 • S * 

2 X « • 1 ERROR «,!> E |2.5,» »*J 

>6 l 

62* 


IF (NDEG « EQ « 1 ) QO TO 100 


>63 



EPS ■ C < 3 > * A 1 | ) / ( C ( 2 ) • A 1 2 I I 

fe . 

>64 

64* 


EE ■ EPS*SQKT|R(2 f 2)yCi2)**2 ♦ 

RC3.31/CI3) - 2«*R(2,3)/IC»2|I* 

>64 

65* 


1 C(3) I ♦ ER A 1 A 2 1 

?.f • 


66 * 


WRITE ( 6 j 9 0 ) EPS,EE 

Kiv 

, ■ • .... - . . _ . - •' 



90 

FORMAT f7?fOX7**EFS ■' i EI2.5,2X, 

» (TRR 0 ^^ ** t E f*. 5 « i ) * F 




GO TO 100 


TV 

u V 


END 

jfti 


ENO of 

compilation: no diagnostics. 

: - ^ 1 ■ " | . 

; : fer 1— 


COMPILATION TIME! 0.95? SEC. 
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POLFIT 


This subroutine fits a polynomial 

pcx) - JZ x k x 


I 


K=o 


to a set of N data points C^t ) Qc ) j l — L - IN/ ,by weighted 

forward regression. The criterion used to determine /t, c is 

2 . 

the minimization of the approximate variance of fit 5 : 

1 ZZ "LiK- 


s 2 ' - 


/v-m -i 


L-i 


where 'VJ'; is the weight assigned to the i-th data point (if the 

L- 

2. I / Z 

variance of the i-th data point is <3^- > then 'Uf’£ ^ / <J£ ) . 

The mechanics of the fit are performed by use of calculated ortho- 
gonal polynomials : 


j3(X) = 23 


where 


k-o 


the R satisfy — i. 


and 


FL/ x) =• ^ Pk '* 


cx; 


A/ 


2 ur- P k tXi> Rj «i) = S hj D k 


t -f 


The fit is continued until all terms whose contributions to the variance 


S are significant at the 90% level (F-test) have been added. 
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CALL POLFIT (X, Y, W, N, C, E, R, M, SIG) 

Input : 

X(I) - abscissa 3C' of the I-th data point 

Y(I) - ordinate ^ of the I-th data point 

W(I) - weight of the I-th data point 

V 

N - no. of data points 
Output: 

i-\ 

C(I) - » coefficient of 30- in final polynomial 

E(I) - 95% confidence interval halfwidth of C(I), i. e . , ~tz OJ 

2 

R(I,J) - covariance of C(I) and C(J) multiplied by t , 

2 

i. e . , t cov (I,J) 

M - degree of final fitted polynomial 

SIG - S , standard error of fit where t is the critical value 
of the student-t statistic for a = .05 (Equal-Tails) and 
N-M-l degrees of freedom. 

NOTE: Because of dimensioning^ the program in its given form is limited 

to N < 21 and M < 4. 



OR •SI POLF I T 

CLE OOO COMPILED BY FOR-V S 7 E 6 ON OH JAN 72 At 08;3HJ03. 


STORAGE USED: CODE I 1 ) 000662! OATA IO) 000*23 


SUBROUT I Nt rOLF * T ( X . Y ■ w , N , c , E , R , M , 5 l G ) g 

FITS A POLYNOMIAL OF DEGREE M ( DETERM ! NEDj-B Y ROUTINE) TO A SET OF 
N POINTS I X ( I ) • Y ( 1 ) ) I » 1 s • , « , N o 

“ w ’ l ' IT'TS THE WEI S H f ~ AS . S ' TSN'Eo TO- Po I N T i • Er* G • FOR AS SUM ^0 CoNSTARI 
ABSOLUTE ERROR, W(I) « 1. , ANn FOR ASSUMED CONSTANT PERCENTAGE 

ERROR, WCl) * l • / € Y < I I ♦ Y f 1 


10 1 
10 I 
10 1 

7 * 

8 * 

9* 

C 

C 

C 

10 1 

1 0 * 

C 

101 

n • 

C 

101 

12 * 

C 

101 

13* 

C 


The polynomial is given by pi 

E ( II IS THE 95 PERCENT CONfID 
R ( I i J ) IS THE COVARIANCE BETA 
SQ U A R E OF THE <f R I T I C AL V A L U E 
THE OEGRiEE IS INCREASED UNjlL 
percent level have BEEN inclu 




DATA V 
DO 200 


I l « 

DO 1 00 
B ( I t J » 
BIS * n 
B ( I » I I 
B ( 5.5) 


SUM I » 
SUM 2 » 
SUM3 m 
V Y o 0 
DO 300 
T * WI 


SUM! a 
SUH2 a 
T a T • 
”SUM 3 “= 
V Y ■ V 
All) « 


D I i 
UII) ■ 
B I 2 I 11 


A Y Y « 


( 1 ) /O • DO/ 
I * 1 , H 


° I . 1 I 
0 • D 0 

0 #uio 

1 a DO 

l .DO 


O.DO 

0.00 

O.DO 

.DO 

I ■ I 9 N 

I) 


sum i t 
SUM 2 ♦ T • 

Y < I ) 

SUM 3 ♦ T 

Y <■ T • Y ( I ) 
SUM3/SUM) 


SUM 1 

SUM2/SUMI 
■ - U < i ) 


0*DO 


T 

T*X ( I ) 







|l 5 *4 

S3* 


VY 3 VY “ A < 1 1 »A { 1 J *D ( 1 » 

155 

54* 


1 ■ 1 ' "J. • 

156 

55* 

500 

I a • i' ♦ 1 ■ i 

157 

56 • 


IF U • GT o 5) QO TO 1100 , • / 

161 

57* 


1 1 ° i - i ■' : 

162 

58* 


SUM 1 ^ 0*00 

163 

59* 


SUM 2 » 0 * DO 

164 

60* 


. SUM 3 a 0 « DO 

1 65 

6 1 • 


SUM 4 a Q.DO 

1 66 

62* 


DO 800 J » 1 »N 

171 

63* 


PI a 1.1)0 

172 

64* 


P2 * 0 « Do 

173 

6 5 * 


WORK a X (J) 

174 

66* 


DO 700 K = 1 p 1 I 

1 77 

67* 


f a ‘WORK “0 ( K ) > * P 1 - V(.K3*P2 

200 

68* 


P2 « Pi 

201 

69* 

700 

P 1 ■ T 

203 

70* 


T a W ( J 1 *P 1 

204 

7 1 * 


SUMl = SUM 1 ♦ T»Pt*WORK J 

205 

72* 


SUM2 a SUM2 ♦ T*P2*W0RK. . $ • . 

206 

73* 


SUM 3 a SuM3 T e Pi ! %/■ 

207 

74* 

800 

SUM4 a SUM 8 + T * Y ( J 1 # . j 

D ( I > » SUM3 F' 

2 1 l 

75* 




U ( I 
V < I 1 
A ( 1 > 


WAY' 

Ayr ■ 


» SUMUSUHJ 
» SUM2/D { 1 1 ) 
: SUM9/SUM3 


» A< I >*A( I ) »Dt I 1 
AYY ♦ WAYY 


900 

1000 


IF (1 o £ Q • 2 ) GO TO 1000 


B ( I « l I ■ -Ull 1 )*B( I 1 , 1 ) - V ( I I ) *B ( 
DO 900 J s 2,11 

B ( 1 . j ) ° B n 1 t J-l ) - m I 1 ) »B ( 1 1 , J ) 
IDF » N-l 

SYN ■ DABS ( V Y - AYY > ✓ I D f 
FP ■ WAYY/SYN 



- VUI)*BI 1^2.1 ) u 

u * bi 1 1 » Ji - vj i n*Bi i-2» Ji 


232 
234 
2 35 
2 3 6 
23? 
240 

88* 
89 * 
90* 
9 l * 
92* 
93* 

i too 

IF 1FP ,GT. F90 f I OF T 1 GO To 500 
MM » 1-1 
M a M M * 1 

AYY a AYY • WAYY . .. .... 

1 DF *» N-MM 

SYN ° 0ABS1VY - A Y Y > / 1 

241 

94* 


~00 Uoo K a l ,MM ‘ . ' ~ 

244 

95* 


T “ 0 • 00 

245 

96* 


DO 1500 J B K i M M . ' . ‘ ' ‘ 

2 50 

97* 

1500 

T ■ T ♦ A(jT*B(J,K» 

252 

90* 

1600 

C (K ) » T ■ - 

254 

99* 


C0£P » SyN*T95( IDF) 9*2 



0* 

1 * 
2* 

DO 2 Uoo k * i»mm 

DO 1900 J ® K,MM 
7 • 0 • DO 


■ ■ ; . . ■ * 



3* 

DO 1800 L a J,MM 





4* 1800 

T a T ♦ B(L*K ) *B(L,J)/D<Ll 





5* 

R ( K » J 1 a T «COEF , 






06* 

107* 

108* 

To?'* 


1 900 
2000 


ElK> a SQRT1R1K,K)» 
S 1 G “ SQRTFSYN) 
RETffffR : — ~~~ 



300 1 10* END 

e n d of compilation: no 

' I ' • . ; • . • ' 

DIAGNOSTICS, 


? ‘ 

. * •/ 


COMF1 La t'l'ON tlM'El 

1/0 REQUESTS! 

I/O WORDS tRANsFEREDl 

I.52& Sec* 
66 
6286S 


:M . 
/’s 

* f i 

pi ‘ ■ 

.yl , 










IV-5 

B. EXAMPLES 

Four example calculations are given. 

1. Data of U. Buck, K. A. Kohler, H. Pauly, Z. Phys. 2M, 180 (1971). 
IWGHT = 1 (assumed constant relative error), so the rms residual 
S for the fit is an estimate of the coefficient of variation of 
the data. The value -S = 0.31% is compatible with known experimental 
precision. 

Examples 2-4 are based upon the synthetic data calculated from (IV-1) 

o 

with I, S^, S^, S^, S^ , given by (IV-2), with e = 8 cpe, r^ = 4 A , 
a. = 0.421359. I.e., 



-a. 

and Y- ^ ( l + & ) where & is a normal distributed 

random variable with mean 0 and variance O"’^’ . In each case 

IWGHT = 1, so the calculated value of S should approximate CT 

A 

2. Y = Y , i.e ., 0=0. Calculated 3 = 0.0098%. 

3. O = 0.5%. Calculated -S = .55%. 


4. o = 0.25%. Calculated S = 0.18%. 






TEST DATA --- EPS • 8 .--- rH « Mi NO ERROR ’ 


NO. OF OBSERVATIONS ... 13 

WEIGHTING SWITCH ... 1 

REOUCEO MASS ... .300000*02 


NO. 


ETR 


Index 


f ; § 

U 




2 

3 

4 


S 

8 


,800000+00 
‘•"9 0 0 0 00 + 0*0 


. 1 00000+01 
• 1 1 0000 + 0 1 
TiTOotJo+W 


.960000+0 1 
T"l 21500 + 0 2 
. 150000 + 02 
. 1 8 1500+02 

• 2160W+0T 


. 1 30000+01 


.253500+02 

.294000+02 


.8 390 90+^) 1_ 
• 7 >'8 7 3 4 ♦ o 1 
. 72 5.4 7 9 + 0 1 
. 678422+01 
T 6367 OB+O 1 


*•+»••*♦* 

S 

. «■ 




.5995 8 4*0 » 


jl04|&7+00 

. 82 30 4 S-o] 

.646667*01 

•5509|H-uI 

‘T"4'629~&j5^6T 


• 394 4,| 7”0 1 
.340Ja6-0l 


. 6 4 1 2 7 2 | 0 l 
.66 7 lH*0'l 


,687979+01 

,705014+01 


TTTWA 9 + 0 r 


.730710 + j|0 



• 1 50000 + 0 1 

• 160000 + 0 1 
• 1 70000 + 0 1 

3 3 7 50 0 + 02 
.384000+02 
.433500+02 

. 536633 + 0 1 
. 5Q9803 + Q 1 
.485533*01 

.296^96-01 
. 260 g$ 7-01 
• 2 3 0 fijl - 0 1 

.7486 99*,! 1 
• 755685**0 1 
.761656+01 

1 1 

• 1 8 0000 + 0 1 

Mini irriii 

.463495*01 

"" • 205 7 ill -O - ! 

. 7^6790*01 

12 

.190000+01 

• 54 1 Sqo + 02 

*443412*01 

• I84&|2-0t 

.771232*01 

13 

.200000+01 

•600000+02 

.425049+nl 

. 1 666|7-0i 

,775097*0? 


POLYNOMIAL OF DEGREE 
|RMS residual ... 


4 DETERMINED 
.97530-04 




TERM 


coefficient 


ERROR 


0 

T 


,81434+01 

-.25576+02 


. 10321-01 
T96T3T7+"oO” 


2 

3 

“4“ 


. 13067 + 03 
-.54394+03 
“.11078+04 


•29574*02 

.36361*03 

»T5322*o4 


A 1 *EPS*RM a .13490+02 (ERROR) 

TpT*rm fHoH a 1 *» 0.421559 ... 


. 17097-01 i 


.32000+02 (ERROR * 

.96126*01 1 


• 40556-0 1 T 


EPS*EPS*RM • .25600+03 (ERROR a 

_ EPS~F R 0 H A 42155 T~W7 


.eoooQ+oi (error * .30107+00 \ 

EPS* • 79998 * 0 l (ERROR * , * 16387+01 ) 



TEST data --- g;PS a 8 • --- RM a 4 k *•*> Sl<*MA * 0,005 


NO, OF OBSERVATIONS ... 

weighting switch .»• 



reouceo MASS •( 

» • » 300^00 + 02 



• 

o 

V 

ETR 

Index 


Y 

»•*«) 

i 

,800000+00 

.960000+01 ,826601+qI ol04j>/ + 00 

»«*#**©**0*e*0©4[ 

2 

3 

4 

.900000+00 
• 100000 + 0 1 
.110000401 

•121500+02 
• 1 50000 + 02 
•181500+02 

.780906+01 
'. 72 1 83 1 +0,1 
.674615+01 

• 82 3@4 5 - 0 1 
,666§||-01 
.550919-01 

• 6 6 9 U 6 5 10 1 - 

.684331+01 

.700827+01 

5 

. 1 20000+01 

.216000+02 

ail> 2 8 4 2 9 ♦ o 1 

.462?]p-dl 

.709! 15 + Cg- 

6 

.130000+01 

•253500+02 

. .598548+01 

• 3944;|.7-ol 

.72 9362+01 

7 

• 1 40000 + 0 1 

*294000+02 

.568245+01 

* 340 106-0 1 

.74304^+01 

8 | 

•150000+01 

• 337500+132 

*531429+01 

•296296-01 

.740894+01 

9 > 

.160000+01 

•384000+02 

• 5 1 1698 + 0 1 

.260417-01 

.758717+01 

10 ^ 

.170000+01 

•433500+02 

.487927+01 

•230681-01 

.765726+01 

11 f 

. 1 800^0 + 01 

•486000+32 

,464704+01 

•20576 1-0 1 

,76896 7~+Ul 

12 

. 1 90000+01 

•541500+02 

.446258+01 

• 184672-01 

.776641+01 

13 

.200000+01 

•600000+02 

.426992+0 1 

i 166657-01 

.778985+01 


RMS RESIDUAL .•• 


•54587-02 




term coefficient error 

>"* • «• * • • * + • a «* • * « • • « • • •"* * « * a » e « « « o 9 a 0 9 


.84388+01 


.44756+03 

-•20993404 


•23523+00 


• 1 &733 + Q2 
•29728+03 
. 16483 + 04 



A 1 + E PS • RM ■ . 13979 + 02 lERROR* .38967400 > 

EPS + RM FROM A! ■ 0*421559 ... •33l6 1 4 02 URROR '<* <>92436 + Q0 ) 


EPS*EPS*RM > •49231+03 (ERROR * ,15748403 1 

EPS FROM A 1 ■0*421559 •«• * 13330 + 02 (ERROR * *43880 + 01 ) 


EPS • *15858+02 (ERROR • 


•Sl9s?+0i 1 



TEST DATA EPS • 8. --- RH ■ 8 , S I GMA » 0*0025 


NO, OF OBSERVATIONS ... 

weighting switch ... 



+ 0! 
, 1 90000 + 0 1 
. 200000+01 


,486000+ 
.5*11500 + 02 
o 60000Q+0 2 


AH 2 57+0l 
• 883228 + 0 1 
.*125656 + 01 


•205761- 

♦ 188672-01 

• 1 66667-0 1 


.768163+0 
a 7 70883 + 0 1 
.776312+01 


POLYNOMIAL OF DEGREE 


RMS RESIDUAL ... 


term 


COEFFICIENT 


3 DETERMINED 


• 18480-02 

ERROR 

TiT*TiTi * i Ti # * * « 


,81921+01 


-.286 I * + U2 
. 18464 + 03 
-• 71388,03 


• 7?R2?-0 1 


« 5 J 1 »6 + 0 I 
* 1 0056+03 
» 55802 + Q3 


Al*EPS+RM « . 13571+02 I ERROR* *l3lS8*00 ) 

EPS + HM FROM A 1 ■ 0*82155? *32192 + 02 (ERROR * 


*31212600 1 


EPS+EPS+RM a *28686+03 (ERROR • .53236+02 1 

EPS FROM A 1 *0*821559 .«* »8S986*0l (ERROR B ' *15719 + 01 } 

EPS • d 10101+02 (ERROR *» *^7869 + 01 1 










V-l 


V. XTREMA 


This main program produces a fit of the experimental variable 


Y = (N - 7 ") v to the model 
8 N 


Y- 


1 


+ 


(A, + 


z 3 


E E E 2 - E 3 

where, after regression the following ansatz is made: 

/c ( - la, er m / j ,c z - 2. A rk j A^-e 


(V-l) 


(V-2) 


which lead to 


a , e r m - 


t 1 - r 


m 


'TT_1* 

2. 

_ /rrjK 
a. 

■ /^3 


y 0 » 

^2. 


(V-3) 


The details of the fitting procedure are described in V.A (FIT). 

The description of the input and output is the same as in Section 


IV. 



3T~X T'R E M a 

000 COMPILED BY Fo m -V S7EA ON 16 SEP 7j A T 19123125. 
ORAGE USED: CODE ( 1 ) On03l7j daTA(O) ODO^A*) 


REaL INDEX (30) i MU 
"FHTETNFT0N E ( 3 0 ) . Y I 30 ) , W ( 30 , , v ( 3 

DATA 41/0*921559/ 

REaD IN EXPERIMENT TITLE 
REaD ( 5 , 5 ) ( T I TLE ( I ) , j « I , 1 fi f 

EOrMat ( 1 8 A H ) 

WR f Te (6,6) * T I T LE c I ) » I => ! t 1 0 > 

F OR* at ( 1 H 1 , 1 0 x , r'8 A R / ) 

READ jN NO. of data POINTS NOB» 


0) .TITLE ( 18) 


Error ANO W c f 1 * 1 • IF 1, RELATIVE 
'and reduced 'mass” Hu. 

NOtE ... IF V IN UNITS Op *M/SEC AND 
of Energy w|ll be i»e-H Ergs, (note 


I.E-27 ERG-SEC) 

RE aD (5,10) nob, I W g H T • M U 
format ( 2 I 5 , F 1 0 . 9 ) 


WR.j Te (A, 20) nOB,1wGHT,MU 

FORMAT I 1 H 0 , ! 0 X , * No • OF OBSERVATIONS 

H I 1 0/ 1 l X «’ REDUCED MASS •••• .EI5 


IES AND EXTRE 
OB) 

® 1 tNOB ) 


WEIGHT SWITCH |iWGHT (IF 0 ABSOLUTE 


RELATIVE ERROR AND;! W ( 1) » |/Y||)t«2|, 


mu I'm i .^2R gram, then uni 
that consShas units of acTI 


»i lO/^lX, * WEIGHTING SWI tic 


1 1H*)/) 

DO 500 I « l t N0B 
E ( I ) > 0 . 5 * Mu* V ( I ) • V ( 1 ) 


wrjT E (&,so» i , v ( i ) ,e ( i ) , Index c i i , v i i i r 

format ( 3X , I 5 » RE 1 5 .6 ) 

S I G ^_F I TJ E , Y , w ,NoB , A l ERM t ERAER ,E2R , ERE 2R. EP S .ER EpS ,RM l ERRM ) 

WR i T E ( 6, 60 )’ S I G" ~ ' ~ 

FORMAT (//JOx. *rms residual OF FIT ••.*,E1S,5/) 

E P S R M « AIeRm/AI 


EE»fRAER/A! 

WRjTe ( 6 , 80 ) AIERM.ERAER, EPSRM.EE 

FOrMat ( / I OX , » A 1 *EpS*RM a * , E 1 2 . 5 , 2X • • ( E RROR ■ « * - ‘ '* 

From a 1 '• 0»R2 IS&9 . • , '• ,E 1 2*5 ,2X , • (ERR 
2 a E2R/EPSRM 

“ rPS?«sQRT( ( ER E 2R /E 2 R ) * • 2 ♦ (EE 


85) E2R, 



12. S»' )»/|OX,*EPS* 


2.5 


WRiTe (6,95) RM.ERRM 

FOrMAT I/lOX,*RM ■ • ,E l 2.5 ,2X , • ( ERROR ■•»E|2»5,*)» » 
GO TO 100 


I 














S3* fc'Ki'n " • ; 

end of Compilation no Diagnostics. 


COMPILATION TIKES 0.731 S£Ci 


I/O 

"RIqufStsI 

51 

}?}••■ # 

I/O 

WORDS TRANSFEREE 

55235 

i 

fS 



V-2 


A. FIT 


This subprogram performs the mechanics of the fit of (V-l) , given 
data E(I) , Y(I), W(I), for 1=1, . .., N , where W(I) is the weight 
of the I-th data point. Initial estimates of and c^ (80 and 5 , 

resp.) are iteratively refined by bilinear regression of the first-order 
expansion 


g - r ) 


/3Y \ 
+ (— ) 

d/Ci 


X 


(£ $ -/>3 ) 


o o 


(V-4) 


where 


9X 

3Y 

9/^3 


( A 4- Az. ^ 4- + 

^ A ' -e" £f 3 


n) 


- 


(A, 




2.A 3 /V ^ 


) 


(V-5) 


The iteration is continued until the rms residual converges to 1%. 


Z - FIT(E, Y, W, N, AIERM, ERAEM, E2R, ERE2R, EPS, EREPS, RM, ERRM) 

returns Z = S the rms residual and 

2 

A1ERMS = a . er , E2R =£r,EPS=e, and 
1 m m 

ERAEM, ERE2R, ERRM, the resp. 95% confidence interval halfwidths of 


these quantities. 



is i r i t • ' 

OoO COMPILED BY Fn^-V S7E6 ON 14 sEp 71 A T 19123127. 

V|jRa 6E UStD: CODE (11 00030*4 1 nATA(O) 000i&7 


l * 


FUNCTION FIT < E . Y • W , N , A l E R M ♦ E R A E R . E 2 R . 

E-R.E 2 R » E P S . r R E P S » R M » E R R M ) 

2 • 

C 

performs iterative fitting of extrema values ... 

3 • 

C 

e ( 1 ) • translational energy of point j 

.1 

R * 

C 

Y ( I 1 * (N- 0 . 37 S)»v< H EXTREMA 

ordinate 

5* 


DIMENSION E( 1 ) »Y(1 ) *VV( 1 1 ,X(30,2) ,C(3) ,FR(R) , AiR) 

6 • 


data A/- 0 . I 65 S 1 D.I 0 S 7 .- 0 . 05 RR, 0,0139/ iCONS/l .6565R701/ 1 

7* 

C 

NOTE - CONS a PI*HcR0SS/2. 

m f 

8* 


EPS a S. 

te 1 " 

9* 


E 2 R a 80. 


10* 


S 1 Cl P ■ 0 . 

Ji. 

1 1 * 


K 0 u N r 0 0 


12 * 

100 

KOljNr ■ KOUNT ♦ 1 

- gfe 

& 1? ' 

13* 


IF (KOUNT ,GT, 10) Go TO RqO 

f . sC . 

Kj£ ’ s? 

«■ '-H 

I R* 


DO 2 Of) I a 1 ,n 

K 

IS* 


X«i*l) e (All) ♦ EpS * ( A ( 2 ) * EPS * ( A ( 3 ) 

♦EPS|§A(R)/E( III /E Cl))/ 

16* 


1 E< I ) 1 /E( I ) 

. 1 
* 

' 17* 

200 

X ( I » 2 1 a E2R *(A(2) +EPS *(2.*A(3) ♦ 3.»EPS • A ER 1 /E ( t 1 ) /E ( I ) ) / 

18* 


t (E( I )*E( I ) ) 

if 1 ■ 

19* 


SlG a BILInR (XiY,W»N,C,eR.RH0,30) 

45- * 

$3. B 

rri -j 

20* 


EPS “ EPS ♦ C(3) ...... 

W f 

fa- 

21 * 


E2R » C ( 2 ) ; 

Hf 

22 * 


DEl a SIG i SIGP 

fi 

2 3* v 


Slf,P • SIG' 


2 R * f 


IF (ABS(DEL) ,GT. 0 • □ t *s 1 G > GO TO 1 00 

i 

2 S • 

300 

AlpRM ■ CONSaC ( 1> 


26* 


ERAER « CONSaER ( 1 1 


27* 


E2R a C0NSaE2R 


28* 


ERE2R a C0NS*ER.(2) 

‘ • 

29* 


ERE p S « E R ( 3 ) 

t; 

30* 


RM » e2R/ C EPSaEPS ) 

k 

31 • 


ERRM . RM*SQRT( R«a(EREPS/EPS)**2 ♦ ( FRE2R/E2R ) a *2 - R.*RHO*EREPS 

32* 


1*ERE2R/(EPS*E2R) ) 

si ; : 

33* 


Fir ■ sig 

4 

3 R • 


RETURN 


IB* 

rOO 

w M 1 T e U,loi ' ” w-*-* 


3A* 

10 

FORMAT t ✓ ✓ 1 0 x • ♦••***FA ILU*E TO CONVERGE 

IN 10 ITERATIONS*/) 

37* 


GO To 300 


38* 


en d 



end Of COMPILATION} 


NO DIAGNOSTICS. 


COMPILATION TIME} 

1/0 Requests i 

I/O WORDS TrANSFeREDI 


o.&rr sec* 

50 
SR 9 7 3 




V-3 


B. BILINR 


Performs a weighted bilinear regression of the model 

Y i ~' /C s+/ c 2,^ti *+ X it 

Z = BILINR (X, Y, W, N , C, E, RHO, NDIM) 

Input: 

X(I,1) - X^, I-th value of first variable 
X(I,2) - I-th value of second variable 

Y (I) , W(I) - ordinate and weight of I-th data pt. 

N - No. of data 

NDIM-row dimension of X 
Output : 

C(J) - c_, , J-th parameter 

E(J) - 95% confidence interval halfwidth of C(J) 


RHO - r 2 g correlation coefficient between c 2 and 
Z = 5 rms residual of fit, calculated from 

9 I ji a ,2 

S = V. ur • ( Y -Y. ) 

to-2 i L h 




S3* Slf,2 a 0,00 

5 *4 * S U M Y e 0 . D 0 r 

55* __ DO 2qo 1 *' 1,N J__ 

56* ‘ T b Y (I ) " C ( I V “ C I 2 ) *X I I ) - C * 3 » • X f !*NOjM) 

57* S I G 2 a S ! G 2 + W(!)*T*T 

58* 200 SUMY a SUMY W ( I 1 • T 

59* C ( 1 > a C ( I I ♦ SUM Y /SUM W 

60* i I 0 r a N**3 

61* S I G 2 a SIG2/IDF 

62* B!|_lMR ■" DSQRT(SIG2) 

63* T a T 9 5 4 IDF1 *b!LINr 

69* E ( 2 > ■ T*DSQRT « A 11 | 

65* “ E ( 3 ) a T * D 5 Q R T < A 2 2 ) 

66 * E ( ! ) a T*DSQRT< <1.00 * 1 A j 1 • SUM X 1 • SUM X 1 «- 2 • D 0 * A 2 i • S UM X l • SUM X 2 * 

67* I A22*SUMX2»SUM* ? » /SUMW 1 /SU M W ) 

6 8* R H o » A 2 1 / D S 0 R T < A 1 1 « A 2 2 1 

69* RETURN 

70* END 


end of co mp il at i on j n o diagnostics. 


COMPILATION TIME: 1.09r 5 £C « 

i/o requests: se 

I To WORDS TrANSFeRedT '58883 




V-4 

C. EXAMPLES 

The. four examples are the same as described in IV-B. 

1. Note that -S = 0.31% 

2. 3 = 0.12xl0~ 5 % (true 0%) 

3. S = 0.66% (true 0.5%) 

4. S = 0.21% (true 0.25%) 



« a- in I'D 03 io* o 


D*fA OF 1 UTBucK, K,A, KOHLeR, *NQ H* PAULY » Z. PHYSIK 2 *1 *♦ 1 00 11971) 


no. op observations ".V. i 3 

WEIGHTING SWITCH ... 1 

REOUcEO NASS ... .392099+02 


NO. 


V 

. ETR 

Index 

y k 



•239000+0 1 

.937699+02 

• 500q00 + 0 1 

-t\ 

.10 6 2 35 + 0 7 


£ 

.207500 +01 

.737336+02 

• 550000 + 0 1 

.106390+0? 



• 108000+0 1 

. 605 2 6 0 + 02 

.600000+01 

. 1 05750+02 



• 1 70500 + 0 1 

.997827+02 

. 650000+01 

. 1 0 0 0 3 IfcjjtO 2 


\ 

• 1 5 9500 + Q l 

.908777+02 

•700000+01 

• 1 02 35^0 2 



.191300+01 

.30191 i+02 

.750000+01 

. 1 0067^102 


it 

. 1 30500 + 01 

.291602+02 

* SOOqOO + O 1 

,99506 ^ 1 0 1 


p 

• 1 20700 + 0 1 

.299985+02 

.850000+01 

.9806 8 Jit 0 1 


H . 

• 1 1 2500 + 01 

.216738+02 

. 900Q00+0 J 

.9703 iff 0 1 



. 1 00500 + 01 

*187009+02 

*950000+01 

• 9 5 3 5 6 2^0 1 



.9^,8000 + 00 

* 1 60065 + 02 

• 100Q00 + 02 

.931700101, 




TT 

13 


.055000+00 


l *♦ 1 B 1 2 + 02 
l 25 1 88 + 02 


» 1 OSoOO+02 
* 1 10000 + 02 


.921375+01 
• 9 0 8 *♦ 3 7 1 0 1 


rms residual 0 p pit ... 


. 30755-02 




If 

4 


a1*EPS*PM ■ .16798+02 (ERROR- .13757+00 » 

ePS»Rm PROM A I ■ 0.921559 . * * .R007r+Q2 (ERR OR ■ .32633+ qO 1 

EPS*EpS»RM ■ .5^1 06+03 (ERROR ■ ,57075+02 ) > 

EPS FROM A 1 »p. 921559 ... .12&I5 +q2 {LOWER BOUND On|eRrOR « .12667*011 

EPS • . 18831+02 (ERROR ■ .2 *♦ I *♦ 2 ♦ 0 1 1 | 

RM- ■ / .15822 + 01 (ERROR ■ ’ ’ A fl * 7 ♦dO’f 


1 ■ ; ! ■ 1 ■ : ' 


• t E 5T -DATA --- NEW Terms - 

- NO ERROR -- EPS 

“: 8 » -- RM ■ 9, 



NO. of OBSERVATIONS ... * 13 ; 

weighting switch ♦ .. i f 

REDUCED MASS ... . 300000 + 02 







,t } 
; 1 ? 

f- 

NO. 

• • • 

V 

ETR 

Index 

4 


1 

2 

. Raooaa+oa 
.900000+00 

• 960000*0 1 
• 121 & 00+02 

. 0390 1 2 + Oj 
•779061+01 

- ~~f- ■ 

.691899+01 
.66 7 9C1&S+0 , 


■ 3 

9 

V 

. IqOOoO + OI 
• li oooo + ot 
♦ 120000+01 

•150000+02 
♦ l 8 1 500*02 
. 216000+02 

•725690+01 
• 678^05 + 0 1 
•636759+01 

.688 1 9.0-+ 0.1 
,7051 0 %+0 1 
.7191 qg * 0 1 


6 

7 

8 

. 1 30000 + 0 1 

; oMOOoo + oi 
•150000+01 

» 2 53500+02 
.299000*02 
« 337500*02 

•599611+01 
• 56M23 + 0 1 
*536693+0 1 

.7 307*ii|+Ql 
.7909f|*0l 
.7987 1%+01 

Jt 

~ 9 
JO 
H 

• 1 60000 + 01 

• 1 70000 + 0 1 

• 1 80000 + 0 1 

• 389000*02 

• 833500*02 

• 986000 + 02 

• 5 0 9 0 1 0 ♦ 0 l 

• 985537 + 01; 
• 963997+0J 

.75569^*01 
• 7 6 1 6 6 ^9 0 1 
. 7 6 6 7 9 0 1 

. i 

12 

13 

. I 90000 + 0 l 
. 200000 + 0 | 

.591500*02 

. 1600000*02 

•993019+01 

•925050+01 

.77123 M^O) 
.7751 o©*0 1 



"T 

rms residual nF 

fit ... 

• 1 22^7-07 




A1+EPs*RM « • 

EPS+RmFROM a [ 

13990+02 (ERROR- .59790-06 > 

« 0.921559 .32000*02 (ERROR - /• 1 9 l 7 1 -n5 

) 


EPS+EpS+RM > 

EPS FROM A 1 ® o • 

♦25600+03 (ERROR * .21102-03 

921559 ... .80000+Ql (LOWER 

) 1 

BOUND ON ERROR « 

.66039-05 


. I 

EPS ■ .80000*01 (ERROR ■ «IS783*oR 1 


V R M ■ ’.‘jOOOo + Ol (ERROR * ’ •12538*08) 




6 


■ ■ 

1 • BOOOOO + UO .960000*01 .82*601 +01 .63128 f+o 1 

2 .900000+00 * j 2 1 5 0 0 + 02 _ .780906 + 01 _ . 6 & 9 0 6 £>§ 0 1 

3 .100000 + 01 . 150000 + 02 " ' . 72 1 03 1 +0 C . 6 8 H 3 3 1 + 0 j 

9 j . 110000 + 01 .t8l500 + 02 .679615+01 .700827+01 

5 . 1 20000 + 01 • 2 1 6000+02 .620929 + 01 .709 1 1 5 + 0 l 

6 . 1 30000 + 6 1 . 2 53 5 00 * 02 .590598 + 01 . 729 3 6 2 ♦ 0 1 

7 .190000+01 .299000*02 .568295+01 .793099|oi 

8 . 150000 + 0 1 .3 37500 + 02_ . 5 3 M 2 9 ♦ 0 | . 7.9 0 8 9 9 ? 0 1 

9 . 1 60000 + 0 1 . 3 8 9 0 0 b + 02 . 51 I 698 + 01 .7587 1 Vlf O 1 

10 .170000+01 .933500*02 .907927+01 • .7&572&|oi 

1] . 1 8000 0*0 1 .986000*02 *969709+01 .7 6896 

TT~ « 1 90000*0 1 .591500*02 .99&258+0J .77*69 

l 3 .200000+01 .600000*02 .92*992+01 .77890 


RMS RESIDUAL oF FIT ... .*6103-02 


Al*EPS*RM » • 1 3*53+02 (ERROR- .32379*00 1 

EPS+Rm FROM Al * 0.921559 .32386*02 (ERROR « j760O8+oO » 

"• " ■ * • ■ ] 

EPS*EPS*RM ■ .30311+03 (ERROR • ,120*3*03 ) 

EPS FROM Al ■q. 921559 ... « ?359 1+Q 1 t LOWER BOUND On ERROR ■ .37312+011 

EPS * .98608*01 (ERROR ■ *79*91*01 1 

R H v '“Tri i7 S *01 (ERROR * 08l82*oT T“ *“ ~ ~~ 





TE STRATA --- EPS . 8. ---RM ■ 0. SlGMA ■ 0.0 025 


NO* OF OBSER V AT i ON^~~i # , ■ ~ “ l3i * ■ ” f 

WEIGHTING switch I [ : / 

REOUcEO MASS ... *300000*02 


NO* 

• • • 

V 

ETR 

Index 

Y || 

• + .««••«• f tv 


1 

2 

.000000+00 

*9o0000*00 

♦960000*01 

• 1 2 1 500*02 

• 038536 + 0 1 

• 7 8 0 3 I 0 ♦ o l 

.60002 8§p t 
• 6^852 l 

1 

. to* 

$■? 

, • !•'£ 

3 

0 

5 

• 1 o o o 0 o ♦ 0 1 

• 1 1 0000 + 0 1 
• 1 20000 + 0 1 

. i 50000*02 • 
. 101500*02 
•216000*02 

• 727 1 00 + 0 1 
•680110+01 
•636197+01 

. 68 9 600 *fjl 
.706871 +Q1 
. 7 1 8 0 3 6 +;| 1 


& 

7 

8 

. 1 30000+01 
, HOOOU + O 1 
« 150000 + 0 1 

» 25 3500*02 
.2*0000*02 
*3i7500*02 

•698302+01 

•566901+01 

•535622*01 

• 7 2 9 l 0 7 ^'Jj 1 
•701217*01 

.7071 83 | 0 1 

• 1' 
SA 

■ t 

9 

lo 
1 1 

i • 1 60000+0 I 

: .170000+01 
. 1 B 0000+0 1 

• 300000*02 

• 933500*02 
» 006000*02 

•511610+0! 
•085033+01 
•060257*01 ’ 

•758589*91 

.762166+01 

.768163*01 


1 2 
13 

. 190000+01 
• 200000 + 0 1 

•501500*02 
• 600000*02 

•003228+01 

•025656+0! 

.770883+Pl 
.77631 2*| | 

? r , ‘ 



RMS RESIDOAL n 

F F I T * * * 

,20500-02 

, li i 

$ ■ 
w 

. i 


ai*eps*Rm ■ 
EPS*Rm FROM a 1 

m 

•13095+02 (ERROR* .98363-01 ) 

a O.02J559 ••• .32013*02 (ERROR a .23333+00 

) 


EPS*EpS*RM » .25010+03 (ERROR ■ .30287*02 

EPS FROM A 1 «*o • 0 2 1 559 *79372 + qI (LOWER 

i 

BOUND On ERROR a 

. 10726 + 01 ) 


EPS '■ .77712*01 (ERROR • 

•25066*01 1 

i 

} 

- 


RM a * 9 20 75*0 1 (ERROR * 

• fl’fOfoTi 





VI- 1 


VI. QFIT 


This main program directs a weighted nonlinear regression of the 


model 

QCir) = V, + Y^[i + H ^] ^ 


(VI-1) 


where 




VE 



+ 


£-' 2 » / (VI-2) 


and after regression the ansatz = £ and (1-8) is made, leading 

to c 4 = * ( YJ 8.o8S) Vi 

^ V H / A-Y/ 

a, = Y 3 Y S /Y, 

% = /^(ir^lzf lz 

e - 


(VI-3) 


The regression is automatically repeated with initial values for 
Y.-Vs- being the values obtained from the first regression, and a new 
adjustable parameter V/ being added to the model: 


Q Cxr) - Vj-U" ^ -h + 


(VI-4) 



VI-2 


After regression the ansatz 

y 2_7T 

» 3 (yi-5) 

is made, leading to 



The required input is of the following form: 

CARD 1: Title information (18A4) 

CARD 2: N, MU, EPS, RM, Al, GO, C6, IWGHT (I5.5F10. 4,E15. 6,15) 

CARD 3 (or more): V(I), 1=1, .... N (8F10.4) 

CARD 4 (or more): Q(I), 1=1, ..., N (8F10.4) 

where 

N - Number of velocities at which cross sections are given 
MU - y , reduced mass of the system 

EPS, RM, Al, GO, C6 - initial estimates for e , r , a. , g , C r 

m 1 o 6 

IWGHT - if 0 data are assumed to have constant absolute error; if 1, 
constant relative error 
V(I) - the I-th velocity observation 
Q(I) - the I-th cross section observation. 

If the scale of the Q’s is not known, C,, C_ and g will be found 

o o o 

finally to contain the same factor, but no other parameters will be 


affected. 


VI-3 


Output of the program includes the values of * j >•••>*& 
with error intervals (denoted parameters 1-6), and the correlation 
matrix of the errors of V, Y 6 ■ 



I — QF-I-T- — — — 

ffD DO COMPILED 


BY FOR-V S7E6 ON 

JlRAGE USEDJ-CODE! 1 > 0002 72 7 DAT A (0 ) -00 734 0) 


1*1 SEP 71 AT 09:22:29. ; 

BL-ANK--COMMON< 21- 000620- 


1* 

— 2 »- 

3* 

4* 

5 * J 

6* 

7* 

— 8 *- 

9 • 

10 * 
11 * - 
12 * 
13* 
-!«#♦— 
15* 
16* 
17 * - 
18* 
19* 
-2 0 *- 
21 ♦ 
22* 
23 *- 
2*1* 
25* 
~26V— 
27* 
23* 
29 *- 
30* 
31* 
-32*- 
33* 
34* 
35 * ~ 
36* 
37 * 
-38*- 
39 * 
40* 
41*- 
42* 
43* 
44 *- 


C 

-e 


A CONSISTENT SET OF UNITS .FOR THIS PROGRAM IS 


•iff 


'it.- 

i,,:- 


100 
jlO - 

i 

120 


C 

c 

c 

c 

30 


' 1 i-E-24-GM--F0R--M AS-S'f~-l-it--14-eR G J S^01^'NeROY^--i^STiWMS— FOP-t«WJT+« 
AND KM/SEC FOR VELOCITY. 

REAL MU . 

DOUBLE PRECISION R< 36 ) »X . P AR~»E R»F<~200 T.F 2 < 200 >Mj§UP<6 lrBtOW<6 I 
DIMENSION XH200>.Q<200).U<200I .PARI6) .ER<6I jpt 
COMMON V < 2001 *E < 200 1 ip©' 

DAT A-BLOH/6*O-*OC-/*0UPy-l-4D6T+^O2T-5Teir»2-rO4-rt»O4V^&6y — 

EXTERNAL XSEfT.DXSECT H 

READ ( 5*10) jt Q (I )» 1-1.19) 

FORMAT CieA47 . - - - : - 

WRITE <6*20) <Q<1)» 1-1 .18) 

FORMAT (1H1 t 10X* ISA 4 ) 

READ <5.30) N »MU » EP S-.R M. A 1tG0-»C6-» IWGHT T" "" 

IF IWGHT r 0. ASSUMED CONSTANT ERROR. IF. I WGHT - 1. ASSUMED 
CONSTANT RELATIVE ERROR. 


1 


fit 


35 


m 


i 


FOR A LENNARD-JONES 12-6 P OTEN TTAL — A 1— =~0 s421 559 »-G0-=- 0 * 1 8 6299 — 

CRUDE INITIAL ESTIMATES FOR EPS AND RM MIGHT BE| 5. AND 4.. RESP. 
FORMAT (I5.5F10.4 .E15.G .15 I 'if 

-W R ITE— <6*3 5) N.MU.EPS.R Rr-AlvG Or e&rTWGHT : 

FORMAT (1 HO » 10 X » * NO OF DATA PO INTS •» 18/1 IX » ’REDUCED MA SS • t El 5 . 6 / 

1 11X. ’INITIAL EPS*. E15.6/11X. ’INITIAL RM’ »El5.pllX .’ INITIAL Al’. 

2 El 5 i 6/ 1 1 X ♦ * INIT IAL- GO * . E15.6/11X . ’ INITTAL €6 ®|e15 .6 /11X . ’ WEIGH TIN 

i 3G SWITCH ’.18/) 

CALL TRIAL ( PAR .C6. Al .GO »E PS .RM1 

REA 0 <5.40) tVtr)« I - l * N 1 

READ <5.401 <Q(I).I=l.NI 
,40 FORMAT <8F10.4) 

J — -DO 200 I - 1 .N — ~~ — * •- — -~ 

W<I) =1. 

IF < IWGHT .EQ. 11 W <I) = 1 ./ <Q <1 1 *Q < 1 1 1 

- E< I)- - -0 *5 *MU*V < I>*V M ) — 

WRITE <6.60) 

FORMAT <//’0 FIT WITH 5 PARAMETERS < NO C8) •//) 

-CALL G A S S 5 9 < 1 » X S E C T » D X S E C T . N * W. 9 .5-.PAR .BUP.BLOW .1.0 -3 . 1 . D -3 • 1 . 0 -6 
1. 20.R.ER.F. F2.X ) 

CALL FINAL <5.PAR.ER.RI 

—WRITE — <6 .70) 

FORMAT <//*b FIT WITH 6 PARAMETERS IWITH C8>*//> 

CALL GASS59 < 1 . XSECT. DXSECT * N.W, Q.G .PAR . BUP.BLOW. 1.0-3. 1 . D-3. 1 . 0-6 


200 


60 


70 


1. 10.R.ER.F.F2.XI 

CALL FINAL <6.PAR*ER 
GO TO 100 

-END- — -- 


.R) 


END OF compilation: 


NO DIAGNOSTICS. 


COMPILATION TIME: 
IVOR £0 0 E5TST -- 

I/O WORDS transfered: 


0.664 SEC. 

51 


55129 


* * 


VI- 4 


A„ TRIAL 


This subprogram converts initial estimates for e , r , a, , g 

m 1 o 

and C, into initial estimates for y , ... , y (denoted by 

b lb 

PAR(l) ... PAR (6) . ) 

CALL TRIAL (PAR, C6, Al, GO, EPS, RM) 

Input: 

C6 - long-range C, 

o 

Al - a^ (I-5c) 

GO - g Q (I- 5a) 

EPS - e 
RM - r 

m 

Output: 

PAR(l) = y 1 (1-8) 

PAR(2) = y 2 (1-8) 

PAR(6) = y 6 (VI-8) 


T R ial : — - — 

000 COMPILED DY FOR-V S7E6 ON 14 SEP 71 AT 09! 22 !4 6. 
AGE USED: CODE ( 1 ) 000102.* D A T A ( 01 000035 


SUBROUTINE TRIAL l PAR t CG » 41 * GO »EPS»RMI 

DOUBLE— PRECISION PAR- — — 

DIMENSION PARIl I ' , > 

DATA HC RO SS/ 1. d5'A» 5919 /.PI/ 3.1*1 159 265/ 


5* - PAR! I I 

6* PAP I 2 I 

7 * PARI 3) 

— ~ 8 -* — PAR f 41 

9* PARIS! 

ID* PARI 4 I 

11* r PARI 6! 

12* r RETURN 

13* { END 

END or compilation: 


8 . 0 8 3 * I C 6 / H C R 0 s S-l * * 0 • *»— : — — — - - — ; — » 

I 2 . • P I * R M I *G 0 * S 0 R T I 2 . * P I • R M *H C R OSSflf S ! 

eps. . ■■ , v . v . fjf ■ : 

4 • *RM •EPS/WCROSS — — — — ' — — ITT 

EPS *P AR 14 I • |g* 

AlfPAR^I ■; . , • ■ § 


NO DIAGNOSTICS. 


COMPILATION TIME: 

1/0 requests: 

-I/O W ORDS- TRANSFEH 


0.395 SEC 
49 




B. FINAL 


This subprogram takes the final values of 

their errors, and outputs values and errors for 
2 

a.. er , e r , $ , and C 0 . 

1 m m o 


Y 1 , . . . , Y 6 and 



CALL FINAL (NP, PAR, ERR, R) 
where: 

NP - no. of parameters of fit (5 or 6) 

PAR (I) Y 1 < I < NP 

l — — 

ERR(I). - 95% confidence interval halfwidth for PAR (I) 

R(I , J) - correlation coefficient r^, between PAR(I) and PAR(J), 
stored as an NPxNP matrix. 


V 


T FIN AL - — r - - :j r- — — ----- - - — : — — 

POO COMPILED BY FOR-V S7E6 ON 14 SEP 71 AT 09:22:49. 


(RAGE USE D-J— CODEC 1 ) 00 05 40* 0 AT A (0 1 000177- 


1* 

— 2 * — - 

. 3* 

4* 

- 5 - 

6 * 

7* 

— 8 * — 
9 * 
10* 

11 * - 

12 * 

13* 


SUBROUTINE FINAL I N P * PA R » E RR *R I 

— D 0 U9L E-P RE C I S ION - P AR-* ERR # R <14- — — — -j~ 

DIMENSION PARC 1 1 »ERR I 1 1 if 

DATA HCR OSS/ 1. 0545 91 9/ *P 1/3. 141 5 92 65/ | If 

" NN " NP * -1 — gr -7 ~ r W~ 

. EPS - PARC 3 1 Li if 

EREPS - ERRC 37 Hr ' f . 

— Al~-r— P A R 1 4 l*E P5 VP A R 15 > — ■— — : — 

ERA1 = A1 *SQRT< <ERRI4 I/PAR C4 II **2 ♦ C E RR f 3 1 /PAf C 3 > 1 * *2 ♦ f ERR < 5 I / 

1 PARC 51 1**2 ♦ 2.*RC3+3*NN)*ERRC3I*ERRC4 )/CPAR^)*PARC4) I -2.* 

2 R( 3*4 * NN I *ERR C 3 1 *ERR C5 1/ C PARC 31*PAR C 5 7 1 -—2 « * ||4 + 4 *NN I *ERR <4 I*— — 

3 ERR< 5 1/ (Par (4 > *PAR C 5» I I ' I ; V 

RM - 0.25* PARC 5)*H CROSS / (PAR 13 l*PARC 31 I 

— ER RM“ RM *S QRT1"'4 ’ir*; C~ER R ( 3 1 7P A R ( 3 l"l *“* 2“ ♦~~("E‘R R C^’l 7 PARI 5 1 1 * *2 — 

1 R ( 3+4 * NN I * ERR C 3 1 *E RR C 5 I /( PARC 3)*PAR ( 51 I 5 f| '% 

CG = HCROSS* (PAR (1 1/8.0831**2.5 

— E RC 6 - 2. 5 *C 6 * E RR C 11/ PA R C 1 1 — ~ — . - - — 

GO = 8.*PARC2 l*PAR (3 1**3. 5/ CH CROSS* HCROSS * C2.*#I*PARC 51 I **1 .51 | 

:■ ERGO = GO *SQRT C CE RR 1 2 1/PAR 12 1 1 **2 + 12. 25* CERRC 3:^1 /PARC 3 1 1 **2 t +2.25# 
-l-IERRt-51/PfRCS 11**2 — *-7-T*R-t-2+2*NNl^ERR121*E:RR|5F/C PARC 2 1 *PARC 3 > h 

2 -3 . *R C 2* 4 * NN 1 *ER R ( 2 1 * E RR ( 5 1 / C P AR C 2 1 * P AR ( 5 1 I - 10 . 5 * R C 3+ 4 *NN 1 » ERR C 3 ! 

3 * EPRC 51/ (PAR 131 *P AR C 5 1 1 1 

- WRI TE C 6 » 1 0 1 — - - - - 

FORMAT C///10X* ’TERM *» 8X . * V AL UE • * 8X t • 95 PERCENT HALFWIDTHV5X* 

1 50 ( 1H* 1/ I 

FORMAT ( 1 0 X * * CG * t 2E 15 . 6/ 10 X * * GO • *2E15.6/10X* * EPS * *2E15 .6 /10X . 
1 * A1 «» 2E15.6/10X. * RM *»2E15.G1 

-A1ERM = 0.25*HCR0SS *PARC4 1 — j 

ERAER = 0.25 *HCROS S* ERR ( 4 ) f 

WRITE 16*301 AlERMtERAER 

— F OR M frT— C / 8 Xr-*- A 1 *EP- S » R M -*rE- l 2 * -&-t E ~ 1 5 * G 1 —7 — — 

E2P - 0.25*HCR0SS*PAR C5 I 'S'-'V 

ERE2R = 0 • 2 5 * HCRO S S * ER R C 5 I 
WRITE C (> * 3 !? I E 2 R » E R E 2 R . 

FORMAT ( 7 8X * * RM *EP S* *2 • *E12 «6 »E18«6t 
IF CNP .LT. 61 RETURN 


-8ETA--= -lv5 *P AR C&17PT 

ERBETA = 1 • 5*ERR ( 6 1 /P I 
WRITE (6*401 BET A* ERBETA 
-FORM A T ( / 1 0 X r * BET A * » 2E15 .& 1 
C8 = CG* BET A 

ERC8 = C8*SGRT ( 1ERC6/C6 1**2 

-*ERBE J T^7-C-CG-*B FT-ftFl- — — 

WRITE (6*501 C8*ERC8 
FORMAT I /10X » * C8 ••2E15 .GI 

RETURN- ,™ 

END 


CERBETA/BETA1**2 ♦2**Rll + 5*NNI*ERC6 


END- OF CO! 


—HO DXAGNOSTXOS.- 


VI-6 


C. XSECT/DXSECT 

Subprogram which calculates Q(v) given y y by means 

1 D 

of (1-6) (XSECT), or 3Q(v)/3y i , 1 < i < 6 (DXSECT). 

1. CALL XSECT (PAR, F, NOB, NP) 

Input : 

PAR(I) - y 1 < I <_ NP 

NOB - No. of velocities v 

NP - No. of parameters (= 5 or 6) 

Output : 

■ F(I) - Q(v) for I-th velocity 

2. CALL DXSECT (PAR, F, NOB, NP) 

Input: 

PAR(I) , NOB, NP. - as above 
Output : 

F (I , J) - (NOB X NP matrix) 3Q(V)/3 y^ for I-th velocity 



I — X-SI-G-T — 

aOPO COMPILED 


BY FOR-V S7E6 ON 1«* SEP 71 AT 09 :22:56. 

RAGE USED: CODEllI 000420; DATA 10 I 00 011-0*- BL A NK— COMMON <21-000620- 


1* 

2*— 

■ 3 • 

4 * 

5 * ~ 

6 * 

7 * 

-a-*— 

9 * 
10 * 
11 * ~ 
12 * 
13 * 

— 14 * — 
15 * 
16 * 

- 17 * - 
18 * 
19 * 

-—20*— 

21 * 

22 * 

- 23 * : 
24 * 
25 * 

——26-*— 

27 * 

23 * 

29 * 

30 * 

31 * 

— 32 *— 


C 

C 

c 

c 

c 

-c- 


c 

c 


33* 
34* 
-35* 
36* 
37 * 
-—38*- 


SUBROUTINE XSECT T PAR *F *NOB *NP I 
-S UB ROUTINE WHICH CALCULATES — CRT) -55 — S-ECT1 0 N5— 0 -1 
OF THE FOLLOWING PARAMETERS ... 

8. 083 * I C6 /HCROSS )• *0 .4 


PART 1 I 
PAR (2) 
PART 3 I 
PAR 14 I 
PART 51 
PARC6 I 
DOUBLE 


r 2. *PI*RM *G0* SGRT ( 2 . *P I* RM*HCROS S/E PSTlifT 


i;,va 

m 


z EPS 

= 4.*A1*EPS*RM /HCROSS 

r4. * EP S * EP S* RM /H CROSS — 

= 2.*PI*C8/CC6*3. I 
PRECISION P AR * F * A # HI # PI7 5* PHI# Z 

DIMENSION PAR C 1 1 * F C 1 1 • A I 4 ~ ■ : — 

COMMON V (200 ItEf 200) 

DATA A/-0. 165500* 0.1057D0* -0.5 440-1*0.1390-1/ 

DA T A HI / 3.-26-7D0 / *P 1 7 5/2-^3 5519 4 4 -5BO/- 

DO 100 I = 1 ! * N 0 B 
Z = PART 3 l/E ( I) 

PHI : -PI7 5 -*■- ( PAR T 4 ) * f A RT 5) *( Al l! *2* IA ( 2) *Z*( A t 3 I 
. ETIII/VTII 


7Z) 

iijlv 

■Iffe- 

•ft 
• *# 


-s$- 


iiS 

ip . 
;..'5 


~PB~ 


100 FTI1 = PART 1I/V(II**0 .4 

-t 1-1 ♦ — PArRT 6 I — — 


Z*A«4I 

♦ PART 2) *0SQRT (V (I I I * ( l^DO ♦H1*ZI*DSIN 




I 11/ 
(PHI 




■ 

it 


RETURN 

ENTRY OX SECT ( P AR . F. NO B# NPI 
E NT R Y P OlNTr WHICH RET UR NS M A TR IX -OF DERI VA 
RESPECT TO PARAMETERS 
DO 500 I = 1 » NOB 

FT IT— = lvD0/VTTI**t^r4 — 

Z = PA R(3I/E MI- 
SUM = ( A T 1 1 ♦ Z* ( A (2 I ♦ ZMAI3I ♦ Z*AC4II 1 l/Et^Il 
PHI r «PI 7 5 * ( PAR ( 4 ) * PfrR( 5T*SUMI/V«T) 

SPHI = OSIN (PHI I 
C PHI = DC OS ( PHI ) 

ROOT V=- SORT (V T IT I — — 

TERM = 1.0 0 ♦ HI* Z 
PROO = ROOTV*TERM 


WITH 


** 

■ B'j 

r i 


39 » 
40* 
41* 
42* 
43* 


- F ( I * NOB I = PROD * SPHI - — ^ — 

F( I* 2 * NO B) = PAR (2 1 *ROO TV* I TERM* CPHI*PAR( 51 • (Af 2) ♦ Z*(2. 
1 3. *Z*A (4 I I )/(E(II*Vf I! I ♦ HI* SPHI l/E (I I 

- PROD -= -PROD *P ART 2 f *CP H I .- fr , — — 

F (I ♦ 3*N0B I = PROD /V (I I 

IF (NP .EG. 6) F ( I*5*N0B I = 1*00/ 


♦ AC 3) ♦ 


500 — FTI * 4 • N OB I- =- P R 00 • S UM / VT1 F 
RETURN 
END 


end of compilation: 


NO DIAGNOSTICS. 


COMPILATION TIME:, 

i/o requests: 

I/O- WOROS T R AN SFERE D f— 


0.978 SEC. 
59 

-5-7954-"- — - 


VI-7 


D. GASS 59 /MAXLIK 


This subprogram performs general nonlinear regression- of NOB obser- 
vations Y(I), with weights W(I), to the model 

Y = Y (©„•••) 

This program is a reworked version of a program of the same name of the 
Univ. Wis. Computing Center. 

1. CALL GASS 59 (10, YMU, DYMU, NOB, W, Y, NP, TH, BUP , BLOW, EPS1, EPS2 , 
EPS3, MAXIT, D, E, F , R, DELZ) 

INPUT: 10 - I/O parameter. If 0, no printing is done. If 1, only last 

iteration is printed. If 2, all printing, is performed. 

A 

CALL YMU (TH, F, NOB, NP) - stores Y(I) in F(I) 1 < I < NOB 

A 

CALL DYMU (TH, DELZ, NOB, NP) - Stores 3Y(I)/30 in DELZ(I,J) 

1 _< I < NOB, 1 <_ J < NP . 

NOB - No. of observations 

Y (I) - I-th observed value 

NP - n, No. of parameters 8 , ... 8 

1 n 

TH(J) - J-th parameter value 0 . Initially should be set to an 

estimate of 8, . 

J 

BUP(J), BLOW(J) - Bounds set upon the allowable values of TH(J): 
BLOW ( J ) < TH(J) £ BUP(J). 

EPS1 - Stopping criterion for relative change in the .sum of squared 
residuals (SSR) per iteration. 

EPS2 - Stopping criterion for relative change in each parameter per 


iteration. 
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EPS3 - Stopping criterion for relative change in SSR from initial 
SSR. 

MAXIT - Maximum allowed number of iterations (e.g. 20) 

F,R - work vectors of size NOB 
DELZ - work space of size NOB; X NP 
OUTPUT: TH(I) - determined value of 6^ 

D(I,J) - (NP X NP matrix) determined value of correlation between 

0 . and 0 , 
i 3 

E(I) - determined 95% confidence interval halfwidth for 0. 

• l 

2. CALL MAXLIK (10, YMU , DYMU, NOB, W, Y, NP, TH, D, E, F, R, DELZ) 

The iterative fitting of the 0 vector is done between entry points 
GASS59 and MAXLIK; from entry point MAXLIK onwards, the statistics of 
the fit are determined. A call to MAXLIK calculates the statistics for 
the given 0 vector, but does not determine the 0 vector by means of 


a fit. 



IT GAS 
00 C 


Irage 


S 59 

OMPILED 


BY FOR-V S7EG ON 34 SEP 71 AT 09:23102, 


USED: CODE ( 1 ) 003321; OATA(O) 0026ZS| 


33 

34 

35 
35 
37 

-38 

39 

40 

41 

42 

43 

44 

45 

46 

47 
49 
49 

-50 

51 

52 


*BL0W»EPS1*EPS2* 






50 


SUBROUT INE GASS 5 9 1 10 »YMU # DYMU »NOB tW * Y «NP »TH*| 

1 EPS-3 » M A X IT t Of E t Pf Rf-OELZl — - 

IMPLICIT DOUBLE PRECISION (A -M -*0-2 I 
REAL Y»W JS* 

REAL STUONT - . _ — ■ | — 

DIMENSION Ytl I f Wfl)*TH( D #B UP <11 t BLOW ( 1 ) # D f 1 1 f E || I ♦ Q < 20 1 1 A< 400 1 • 

1 P(20ltPHI(2G1 t TB ( 20 ! «HTS<2CI*F<ll»R(ll» DELZt 1 1 j> | 

— D A T A FLAM/ l"i0”2/» FNU / 1 t 01/~ * If- 

N PS 0 = NP*NP , fr 

GA = FLAM & • 

IF <10 .EG. 0» GO TO 50 ' ; • 

WRITE f 6 »l ) NOOfNP 
FORMAT («1 ENTRY TO G AS S 59 V /I OX * ♦ F IT TO* fI8t3Xf? OBSERVATIONS* /10X 

I t ‘THERE- ARE * » I8 # 3 X » ’ P A R A MET ERS XN~TH E~M00£lr*7 ) % — 

WRITE < 6 f 2 1 ' if. t 

FORMAT C // 5X » * IN IT IAL PARAMETER VALUES *//l5X » ®L^JIER BOUND »f 9Xt *PAR§ 

1 A M E T E R V A L UE * » 5 X f • U P P ER - B 0 UN U * 7 5 X~r7iH- **■ 

WRITE ( 6 1 2001 I (I»BLOWCII»THIIl » BUP f II f 1=1 eNP J 
ITER =1 

— 5 S 0-=-G •- DO- — 


■SB 


I-- 


m 


~?T 


.ff. 


l 


CALL YMUf THfFtNOB tNPJ 
DO 100 I = If NOB 

- - — R(I ) r~YI I 1 *-F(IJ — -r— — — j r - 

100 SSQ = SSQ ♦ W(I)*R(II*R(I) ■< 

KOUNT =1 j 

_ _ — . — s s-qIn T —r— SS Q — : j 

SSQB4 = SSQ f 

IF (10 .NEo 01 WRITE (6.3J SSQ ' 

“3 F'ORMAT— f// 5Xt*INI TI A L S UMOFSGUAR ES* f D 20-rl 0 1 

IF (SSQ .LT. 1.0-501 RETURN 
150 CONTINUE 

— - I F (10 vGE J" 'Z ) 'WRIT E'“ 1'6 ? *(T“TT EITtKIWNT 

4 FORMAT C///5Xf * AT ITERATION* fI5tlOXf ’CUMULATIVE NO. FUNCTION CALLS 

1 MADE* 1 1101 

G A =" GA ♦DSQRTCSSQ/SSQBRI /FNU — — — * — — : 

SSQ84 = SSQ 
1NTCNT =0 

C ALL DYMU ! T H • DEt ZtNOOtNPI — " — — 

NI = - NOB 
DO 650 I = ItNP 

•= -NI ♦ -NOB — : — — — — : t — 

= ( I- 1 ) * N P 
= I - NP 


500 


NI 
II 
I J 

N J-~ NOB— 

DO 550 J = If I 
NJ = NJ ♦ NOB 

IJ : IJ 4 NP 

JI = II ♦ J 
SUMD = 0. DO 

-00-50 OK=-lfNOB — — 

SUMD = SUMD ♦ OELZ(NI*K J *DELZ( NJ+K I *W I K I 
D( JI ) = SUMD 




60 * 
SI* 
— 62 *- 
63 * 
64 * 
— 6 5 *- 
66 * 
67 * 
--68 * 
6 9 * 
70 * 
— 71 * 
72 * 
73 * 
-74* 
75 * 
76 * 
— 77 * 
78 * 
79 * 
— 80 ♦ 
81 * 
82 * 
— 83 * 
84 * 
85 * 

-- 86 *- 

87 * 

88 * 

— 89 *- 
90 * 
31 * 
-92 * 
93 * 
94 * 

— 95 -*- 
96 * 
97 * 

— -98 *- 
99 * 
100* 
-101 •- 
102* 
103 * 
104 • 
105 * 
106 * 
- 107 -*-" 
108 * 
109 * 


- DITJ)- x SUMO 

Q (I ) = 0.00 
00 600 K = 1 »NOB 

GUI = G( 1 1 ♦ • W ( K I * DEL Z C NI +K I* R C K ) 

II = II ♦ I 

IF (Dim i.L 7 . 1 . 0-501 Dill) = l.D -50 

-EC -*>-=— 0 S 9 RTH-D IHH 

NI = - NP 
DO 710 I = 1 * NP 

NI - NI ♦ NP— ; — 

IJ = I - NP ! . : 

WTSCI ) = 0 • DO •: 

-0 0 - 710 J -=- 1*1 

IJ = IJ ♦ NR 
JI = NI ■* J 

TEMP = -DMJI 

AC IJ) = TEMP 
A ( JI ) = TEMP 

- C ALL’ S YM QR-t A-fNP^NPrPTr-OE tZ^ 1 

NI = -NP 
DO 730 I = 1 *NP 

NI = NX ♦ NP : ~ 

DO 730 J = |#NP 
JI = NI ♦ J 

WT S(J) = W T S ( j ->— * A 0 S^ OE LZ f J rt-t^Pt lt — 
MASTER =1 
DPROD = 1.00 

ND = NF t - — : 

DO 770 I = l.NP 
NI = NI ♦ NP 

~IF"iWT SC r) ~T 0 T * — WT STMA^S TE ft ) ) M A STE R = 1 
DPROD = DPROD*D(NI*I)**fl./NPI 
IJ = I - NP 

DO 770 J = 1*1 — 

IJ = IJ ♦ NP 
JI = NI ♦ J 

A <~D) =— DCIJ)/CE<-D* E C J)) — - 

AC JI ) = AC IJ) 

IF CDPROO ,LT. 1 . 0 - 50 ) DPROD = 1 . 0-50 

WTSMST = WTSC MASTER) . -- - - 7 ~ — 

DO 780 I = 1 * NP 
II = CI- 1 l*NP ♦ I 

P CD-OCI) /EC I > — — — 

PHI Cl ) = PCI) 

IF CDA 8 S CWTSCI) I .LT. l.D- 50 ) MTSCI) = 

WT S CD = WTSMST/WTS M ) 

ACII) = ACII) ♦ GA*WTSm*OPROO/OIIII 
CALL MATIN C A » N P * P * 1 # DE T *jP I VRA T I 

STEP = 1 . DO • : i <• " ; 'V • 

SUM1 = 0 .00 

SUM 2 = O.DO * 

SUM 3 = 0.00 
DO 790 I = l.NP 

-SUM 1 = PCI ) * P HI CD+-SUM 1 - 

SUM? = SUM 2 ♦ PCI )*PCD 
SUM 3 = SUM 3 ♦ PHIC I) *PHIC I) 




110 * 
111 * 
112 * 
113* 
114* 
115 • 
- 116 *- 
117* 
118* 
119* 
120* 
121 * 
-12.2*- 
123* 
124* 

- 125* 
.126* 

127* 
12 8 * 
.129* 
130* 
131* 
132* 
133* 
-134*- 
135* 
135* 
137* 
138* 
139* 

140*- 

141* 
142* 

- 143* 
144* 
145 * 

—145*- 
147* 
-148* 
149*- 
150* 


790 


i 


= 1.D0 


800 

810 


820 

7 

-83 0 


— 340 


; 8 

18 50 

i86 0 

8 70 


PHI(I) - Fill - 

TEMP = SUM1/0S0RT <SUM2*SUM 3) 

IF (DA8SCTEMP) .ST. l.DO) TEMP 

TEMP z 57.29578D0*DAC0S (TEMP) — 

IF CIO ,LT. 21 GO TO 800 
WRITE (6*61 DETt^IVRATrTEMP 
— FORM A T it 9X » * 0 E'TERft^N'AN-T*-r02tTvtrO-r9'X‘r*'^A J T'fi) 
19X » 'ANGLE IN SCALED COORDINATES* *F 10 • 5 ) 

00 810 I = 1 * NP 

PCI ) = PHIC 1 1 * STE P/E C I) 

ISW z 0 

DO 820 I = l.NP 

— T 8 C I > - THC I) ♦ PC I) 

IF CTBCII .GT. BUPCI) .OR. TBCI) . LT • BLOW t 
IF CIO .LT. 21 GO TO 830 

WRITE (6*7) • — 

FORMAT C /9X * ’TEST PO IN T P AR AMET ER VALUES*) 
WRITE (6*2011) CTBC I) *I=1»NP) 

— IF-MS W - • E Q vfl ) GO- TO- 860 : 

SUMB r 0. DO 

CALL YMUCTB#F»NOB»NP ) 

-KOUNT = KOUNT * 1 — : ‘ 

DO 840 I z 1 * NO 8 
R (I ) - YCI) - PCI) 

-SUM9- z-SUMB- f - W C H *R CT)*RM-)-— 

IF CIO .GT. }.) WRITE (6.8) SUMB 


1 ) ) 


* 'Z 




- i 


F; 




♦ EPS1)*SS0I 
.AND. GA .GT. 


-8 80 


890 

900 


FORMAT C/9X * *TEST POINT 
IF C SUMB • LE<. (tiOO 
IF (TEMP .GT. 3.01 
STEP z 0.5D0*STEP 

"INTCNT "Z— INTCNT *~T — : 

IF (INTCNT .GE. 36) GO TO 
GO TO 800 

GA z GA*FNU- 

INTCNT z INTCNT ♦ 1 
IF (INTCNT .GE ■ 36) GO TO 

-GO-TO— 760 : 

DO 900 I z 1 ,NP 
THCI) z TBCI) 

-SSG-Z SUMB - 

IF CIO .GE. 2) WRITE (6*10) 


SUMOF SQUARES **020.10) 

TO 


GO T0890 
l.D-50) GO 


~w 


880 


980 


980 


GA 


•151* 10 FORMAT C/9X*»AT THE END OF THE REGRESSION* GAMMA IS**D20.10) 

15 2 * — IF C E P S2— yL E~«— C«rOO-) — GO — T 0 920 

153* DO 910 I z l.NP 

154* IF C DABS CPU) ) .GT. OA 8S ( TH ( I ) * EPS2) I GO TO 920 

— 155* 910 CONTINUE : — * — : — 

156* WRITE ( 6*11 ) EPS2 

157* 11 FORMAT (/ *OITERAT ION STOPS . .. RELATIVE CHANGE IN EACH PARAMETER I 

158 * — IrS— tES-S— T-HAN-*-r020 .IrQ-) ; 


159* GO TO 999 

160* 920 IF ( DABS (SUM B -SSQB4) .GT. EPS1*SSQB4) GO TO 930 

-16!* WRITE 16*12) EPS1" : — 

162* 12 FORMAT (/* ITERATION STOPS ... RELATIVE CHANGE IN SUM OF SQUARES I 

163* IS LESS THAN* *020.10 ) 

— 164* GO-TO -9 95 — — — — 1 — * 

X * 01 AG NOSTI C * THE TEST FOR EQUALITY BETWEEN NON-INTEGERS MAY NOT BE MEANINGFUL. 
,165* 930 IF (STEP .NE. 1.00) GO TO 970 


DO 950 IDUB =- l*-5 i — 

DO 935 I = l.NP 

TBII1 r THIII ♦ Pill 

IF ITBII) . GT » BUPfl) *ORs—TBM) 

CONTINUE 

SUMB = 0. DO 

G ALL YMIH TB.DELZ.NOBtNP) — 

KOUNT = KOUNT ♦ 1 
DO 940 I = 1 » NO B 

SUMB =~ SUMB ♦ W(I)*tYM) — ~BEtZ Ml )♦ * 2 
IF ( ( 1 .DO ♦ EPS1 ) * SUMB .GE. SSQ) GO TO 
DO 945 I = l.NP 

T H l I )-TBtII 

DO 950 I = 1 * NOB 
R( I) = Y (I) - DELZIII 

F (I I - DEL ZCI) 

IF (10 • L T • 121 GO TO 9G0 
WRI TE (6*7 I » 


BL 0U I I M 6 0 -10-970 


R( I) = 
F (I I = 
IF (10 
WRITE I 


WRIT E — ( 5 *203 
WRITE (6.81 
SSQ = SUMB 

CONTINUE - 

ITER = ITER 
IF (SSQ .LEi 


( 5 r 2 0 l-l-h-T-T B ( 1 1 -r^ s-tr-NP- 
(6.81 SUMB 


IF (SSQ .LEi . E PS 3*SS QI NT ) GO TO 991 fi 

j f — (17 r R jlAXITIl 50 .150.9 92 $ : 

980 WRITE (6*13» 

13 FORMAT f///*0****SUM OF SQUARES CANNOT BE REDUCED TO SUM OF SQUARE 

— IS AT LAST- I TE R AT I ON-TER MIN AT IfrN 06 C UR S * / ) 1 

GO TO 999 

990 WRI TE (G *1 4 I DE T .FIVRA T 

14 F 0 RM A TI / • O R AT 10" 0 F~ PIVO'IS—IS — TGO—SnAtlrV-*— TTEI-ER M IM^ffT — •nr5X^"O2O*»-i'0T* 

1 5 X * ’RATIO OF PIVOTS* tD20»10) ! 

GO TO 939 

991 WRITE (6. 157 EPS 3 

15 FORMAT ( /’OITERATION STOPS ... CHANGE IN SUM OF SQUARES (RELATIVE 

IT 0 INITIAL VALUE ) IS LESS THAN *>DZO. 10 * % 

GO TO 999 1 : 

992 WRITE ( 6 1 1 G I MAXIT 

16 FORMAT (/» •••♦•ITERATION TERMINATES— NO. OF ITERATIONS EXCEEDS** 

- liioj - ■ ----- •' — 

2001 FORMAT (19. 3020. 10) 

2011 FORMAT (5X.6D20.10) 

2012 — FORMAT ( / 110. 6020 .-10/MOX *602 G-wl-Of* : — 

GO TO 999 

ENTRY MAXLIK f 10 • YMU. OYMU .NOB • W. Y . NP.T H. D » E. F.R . DEL Z I 

— N PS G NP *NP — — — — - — 

999 SUMW = 0 .DO 

CALL YMU( TH.F.NOB »NP) 

SS 0 =- 0. 00 ■ — • 

IT (10 .OT. 1) WRITE 16.1021) 
no iaoo t - l.Non 

- SUMW - S SUMW ♦ W(II - - — 

R( I) = Y(II - F( I) 

IF 110 .LE. 1) GO TO 1000 

TEMP-=OSQRT(W(I))*R(I) 

WRITE 16.1019) I*Y(I) »F (I) *RII I.TEMP 

1000 SSQ = SSQ ♦ W (I l*Rf II*R(I ) 


1000 



200 0 


1015 


2100 


S » 1 01 5 1 
MATRI X’/T 


101 1— FORMAT ( 1 HI *5 X ** NO . ’ *8X • ’0BS£-R-VE-D-’-*-7X*-’-F-INAL—M0DEL— V-ALUE-5 * * -9 X* 
l’RESIDUAL ’*8X*’WEIGHTE0 RESIDUAL* / 5X * 100 1 1H* I / I 
1019 FORMAT <19*4020.101 

— r- - — CALL DYMUC TH.OEL-ZtNOB *NF I 

NI = -NOB ... 

DO 1900 I = 1 *NP 

' — tN i-= ni-^-nob ; — — — 1 

II = II -1 MNP 

NJ = - NOB j& 

r 1 j r I NP - — — ” * - - 

00 1900 J = 1*1 : f. 

NJ = NJ * NOB B 

I J -r f J- ♦ N P 1 M 1 


JI : II ♦ J |f 

SUMO = 0 . DO "I, 

DO 1800 K - 1*N0B~ H * 

1800 SUMD = SUMO ♦ W<K)*DELZ<NI+K) •DELZ<NJ*K) & ' it 

A (I J ) - SUMO $$ ' ..if 

-1900 — A< JI) --SUMD 1 r *r 

2000 IDF = NOB - f<P fS 

IF <10 .NE. 0) WRITE! 6* 1 01 51 *) 

jl015 ~ F ORMAT < ♦ 0 C 0 R R E LA T 10 NM A T RTX*-/» — — — — “ 

DO 2100 I =■ ltNPSQ 

2100 Dill = A I 1 1 h 

— C ALL - M A T IN -< D *N P* A-r&*0 

DO 2200 I r l.NP if 

II = ( I- II *NP ♦ I , ; • ’ ij& 

- IF ( OCIII' *LT. -lTD-Sd-IHfli = t »0 » S0 1~ 

2200 EIII = OSQRTIOlIin $ 

NI = - NP 5 $ 

_ 0 q— 2 UOO -I r-l-wNP : : 1 — 

NI r NI . NP 

IJ = I - NP i 

DO 2300 J -1>I : 1 — — 

JI = NI ♦ J J 

IJ : IJ ♦ NP 

D < J 1 1-0 < JI 1/ <EM )*E<JM — 

2300 DIIJI = D ( JI I 

2400 IF <10 .NE. 01 WRITE <S»2007) I* < Df NI* JI * J-l * I) 

2007 FORMAT ~<5Xt ’ROW* *15 *9F12 .6 /« 13 Xr9F“12r6 M ~ “ — — “ 

IF <10 .NE. 01 WRITE IS»31) DET * PIVR AT 

31 FORMAT 1/ * ODE TERM IN ANT* » D2 0. 10 »10X**RATI0 OF PIVOT S* *D20.10 I 

j F-f-l DF-iEQ-.-Ol-G0T0-3OOO r 

SDE V = SSQ/IDF 
T F ACT R = STUDNT < IDF 1 

IF <10 .NE. 01- WRITE- 16 91014 l-SSQ^SDtVffTDF 

1014 FORMAT I /// 5X . * SUM 0 F S3. R ESIDUALS * *020.10/ 5X* ’MEAN SO. RESIDUALS 
1* *D 20. 1 0* 5X * * < FOR * » 18 »3X * * DE GREES OF FREEDOM) * I 

— SO E-V---0 SG R-T-T-SO E-Vl 

IF MO .NE. 01 WRITE 16*20081 SDE V *TFACTR »SUMW 

2008 FORMAT 1 5X » * ST AN DARD ERRO R OF EST IMATE’ *020.10* 10X * ’CRITIC AL STUDE 

- 1NT-T VALUE * *F10 .3/5X * ’SUM OF -WEIGHTS * *020*10? 

IF <10 .NE. 01 WRITE <6*20091 

2009 FORMAT <// /5X. ’NO .** 7X * ’PARAME TER VALUE * *4X * ’ST AND ARD ERR0R’*7X* 


2200 


2300 

2400 

2007 


DO 2500 I r l,NP 
TEMP = EfIMSDEV 



-280 • 
281 * 
2 82* 
283* 
284 * 

285* 

285*- 

287* 
288* 
289* 
290* 
291 • 
-292*- 
293* 
294* 
-295 *' 
2 95* 
297* 
“298* 
299* 


2 50 0 
3000 


3100 


3200 


2003 


— El ■!.)—= r l EMP *7 FAGTR — — — — 

IF (10 .NE. 01 WRITE 16 *20011 I* TH 1 1 ) * TE MP *E C 1 1 
IF (10 .EG* 0 « 0 R • PIVRAT .LT. 1*0-121 RETURN 

DC 3100 I = 1 »NP — ~ r-r~ 

II = (I-1)*NP ♦ I 

oti » = Aim !% 

— G A L L— SYMQR (A* N P *NP * P* DEL ZF - — rt : - 

NI = - NP # 

00 3200 I = l.NP 10 

NI = NI ♦ NP : §§" 

At NI ♦ I I = 9(1) . 

IJ : I - NP 

— 0 0—32 00 — J- — Irl f-r, 

IJ r IJ ♦ NP ft 

A(IJJ = A<N|I*J 5 if 

■ WRITE -16 *20031— ’ r if 

FORMAT I//9X* 'EIGENVALUES OF MOMENT MATRIX-FINAL ANALY' 
WRITE (6*2011) IP(I).I=1*NP! 

— R ETURN : — * ; 

END h 


--E ND- OF- COMPIL ATI ON I 


1-0 1 A GNOSTIC! 


GO MPIVATH)N--T-I14e-lf 

s i/o requests: 5 

I I/O WORDS transfered: 


124 
945 57 
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E. SYMQR 

This program is a reworked version of a Univ. of Wis. Computing 
Center routine of the same name which computes the eigenvalues and eigen- 
vectors of a N X N symmetric matrix A by means of the Q-R method. 

CALL SYMQR (A, N, NDIM, EIG , VEC) 

INPUT: A - matrix of which eigenset is to be found. (Only the lower 

triangle of A is used) 

N - size of A 

NDIM - dimensioned size of A 
OUTPUT: EIG(I) - I-*th largest eigenvalue of A 

VEC (J , I) - (NDIM X N matrix) - J-th component of eigenvector I 



SYMQR 1 

0 COMPILED PY FOR-V S7E6 ON q 6 SEP 71 AT iOj*l31H2, 
AGE USED i CODE ( J ) Q0IS6II D A f A ( 0 | 0 0 1 D 5 2 


^SUBROUTINE S Y M (JR ( A , N . N D I M , E t G . V £ C ) ; 

USES THE QH ALGORITHM ON SYMMETRIC MaTRTX A To OBTAIN EIGENVALUES 
NOTE THAT ONLY LO*E« TRIANGLE OF A IS USED (AND pESTROyCD) 

E1G - VECTOR Of. EIGENVALUES I N DESCEND I NG -OROeR^ - 

vec - matrix of eigenvectors in, order of eigenvalues 

IND - ERROR RETURN INDICATIR 

i f . o. normal, -L. sum.-Of - eigen values- not— trace , 

IF 2 SUM OF SQUARED E|G. NOT EQUAL To NORM ||j ■- 

IF 3 BOTH OF ERRORS ;j?§| ' ft 

IMPLICIT DOUBLE PRECISION |A-h, 0-2) l}£|. 

DIMENSION A(l) .GAMMA (SO) , BE Tat 50) .BETASQ(SO) t> E l « 1 t »WI*49) ,VEC(l ) . 
1 IPCR9) ,Q(H9) , } IpoSV(50) .IVPOSISO) ,10RD(S0) % 

1(1 ) ) . (PC 1 ) .BETA! U » • ( Q ( I I *BEtaI I I ) & 

I NO « 0 

..IF • (N . E « . 0 ) GO TO - 560 — ~~ 

N l o N “ l 


N2 » N 
.ENORM - < 
TRACE * 


- 2 

0*D0- 

1 o»oo 


■ - NDIM 
l ) Q J » I , N 
a NJ ♦ NO I M 
" NJ ♦ J f 
1 00 — I — a--J * N- 


IJ « NJ 
ENORM a 
TRACE a 
ENORM ■ 
ENORM a 
-GAMMA ! I ) 


ENORM ♦ 
TRACE ♦ 
ENORM - 
ENORM ♦ 
•--AM-)- 


A ( I J ) * • 2 

A (Jj> 

0*SOo*A(JJ)**2 

ENORM 


IF ( N2 ) 280,270, 120 
NRN ■ “NO I M 

-DO 260 NR ■ 1 , N 2 

NRN » NRN ♦ NO |M 
. I SUB » NRN ♦ Nr ♦ I 

B ■- A ( I SUB )-— 

5 a 0,00 

DO 130 I a NR ,N2 
S a S ♦ A ( NRN J > 2 ) **2 
A ( I SUB ) » 0,00 
IF (S oLE. 0 • D o ) GO TO 2S0 

_5.-a_S~.4- .B • B_._ 

SGN a 1 .00 

IF (B ,GE. 0*00) GO TO 16 
SGN ■ •! .DO 
SQRTS » DSQRT(S) 

D a 0,5Do*SGN/SQRTS 
_T E M P— a-0 S Q R T-CCT, &DT3-^-B • 0-) 
WINR). a TEMP 
Al ISU b ) • TEMp 
D ■ d/temp 


GO TO 





S3* 

5 4 • 
55 * 
56 * 
57 * 
58 »- 
59 * 
60 * 

6 1 *- 


170 


62 * 


63 * 

180 

.. 64 * 

— 

65 * 


66 * 

190 

67 * 

— 

. 68 * 

200 

69 * 

210 

— 70 * 

220 

. 71 * 


72 * 

230 

73 * 


74 * 


75 * 


-7 6 * 



77 * 


78 * 

2 40 


79 * 
80 * 
61* 
8 2 *~ 
83 * 
84 • 
85 * 
86 * 
87 * 
... 8 8 *- 
89 * 
■ 90 * 
9 |* 
92 * 
• 93 * 
- 9 H*- 
95 * 
96 * 
- 97 * 
98 * 
99 * 
l 00* 
101 * 
102 * 
103 * 
104 * 
105 * 
106 *- 


250 

260 

270 


280 


290 


300 


310 


1+ 2) 


210 


107 * 

108 * 

• 01 AGNOSTIC* 


B a'-SGN*SQRTS 
00 1 70 I « NR . N 2 
TEMP b 0 • A ( N R N ♦ 

W(iM) 8 TEMP 
A ( NRN ♦ 1* 21 • TEMP 

-WTAW * 0.00 — 

00 220 I * NR.N1 
SUM » 0»DQ 

00 180 J ■ NR * J 
1J a j*ND 1M ♦ 1 ♦ 1 
SUM * SUM ♦ A( 1 J ) * tt ( J ) 

-it a -- l + 1 

IF ( N 1 .LT» 11) GO TO 
DO 200 J ■ 1 l * N 1 

1 J * I *N0 I M ♦ j ♦ | — 

SUM ■ SUM ♦ A ( J J)*W( J) 

PM 1 a SUM 

-IN TAW — o W T A IN- ♦- SUM* *1-1-1 

DO 230 I * NR.N1 
Oil) » P 11) “«TAW*VM I ) 

DO 24 0 J » NR.N1 

QJ « <J(J> 

W J « VM J 1 1 

-DO 240-1 a--J.Nl ~~ — — r 

1 J ■ J*NO 1 M * l ♦ | ' 

A 1 I J 1 « A ( 1 J ) - 2*00*(W< I )*0j * NjlQ(ll) 

BET A 1 NR 1 « B ...... .. 

BETASQ(NR) a B • 0 
GAMMA(NR*1) a AUSUB ♦ NO 1 M 1 

I J » n2*NDIM-*- N— ; —— 

B a A ( 1 J ) 

BETA 1 N-l 1 a B 

BETAS « < N - 1 ) a 0*B- 

GAMMA ( N 1 • A 1 ! J ♦ NOlM) 

BETASQ(N) a 0.00 

DO 300-1- ■ ! .N 


m 


$ 






! «t! 


. c .. 


♦ I 


IJ • I - NO i M 
DO 290 J ■ 1 .N 
I J » 1 J ♦ NDIM- 
VEC ( 1 J 1 a 0 » DO 

l J • 1 l-l 1*N01 m 

-VE CM J 1 « t ♦00— 

M a N 

SUM a 0 • DO 

NP A S- a- 1 * 

GO TO 400 

SUM ■ SUM ♦ SHIFT 

-COSA — a — 1«D0 ; : . 

G a GAMMA! 1 1 - SHIFT 
PP ■ G 

PPBS ■ PP*PP * 8ETASQM1 — — 

PPBR a DSQRT 1 PPBS ) 

NJ a - NOlM 

-00 — 3-70— J — *— 4-* M — — 

NJ a NJ ♦ NO I M 
COSAP a COSA 

THE- T ES T- F 0 R- E-0 U A E |TYB ETWE 6N— N ON- 1 N T E GER S- 


M Ay N O-T-BE— ME AN l NGF U L , 


I 


IF (PPBS ♦ N £ • Q, DO) GO TO 3*0 
S 1 NA * O.DO 

SINA2 • 0*00 ' : 

COSA ■ I .DO 
GO TO 3&0 

S l N A B E T A I J ) /.pPUR «■ — i 

SINA2 « BETASQj JJ/PPBS 

COSA a PP/PPBR , 

. NT • J ♦ NPAS — - jj 

IF (NT • GT, N) NT - N j>l 

DO 3*10 * ■ I » NT p\ 

-I J * -N J ♦ 1 -1 : : Ui 

TEMP ■ COSA •VEC(IJ) ♦ S!NA*VECtlJ ♦ NO I M ) 

V£ C ( 1 J ♦ NDIM) s -S I NA* VEC ( I J ) ♦ COSA* VEC ( I J+ NO 

VEC(IJ) = TEMP. _____ 

DIA » GAMMAIJ+jJ - SHIFT 
U « S I NA2 ® ( G 4 0 I A ) 

_gamma i j ) »._6 --p. — u — — j 

G * DIA - u S | 

PP » DIA*COSA - SlNA*COSAP*BETAt J) 

IF (J .NE. M) go To 360— — — 

BETA(J) a S I N A * P P 
BETASQ(J) » SInA2*PP*PP 

GO- TO -3 BO — — — 

PPBS ■ PP«PP ,♦ BETASQ(JM) 

PPBR » DSQRT(jPpaS) 1 

BETA! J) « S I N-A *pPBR - - — 

BETASQ(J) - £lNA2«PPB$ 

GAMMA(MM) » G 

NPAS ■NPAS* — J — 

IF (BETASQIM) .GT* 1 • D - 2D GO TO H 1 Q 
E I G ( M ♦ 1 ) ■ GAMmA<M*1) ♦ SUM 

-BETA (til •-* O.DO f 

BETASQ(M) ■ 0.00 . 

M a M - 1 ' • $ 

. I F- < M - ♦ £Q . - 0 ) -GO- T0-H30- — — — 1 

IF (BETASQIM) , lE . l.D - 21) G'O TO 390 ) 

A 2 • GAMMA (M* | | 

R2 " 0 « 5 D 0 • A 2 

Rl a U.!»DO*GAMmA (Ml 

R 1 2 » Rl ♦ R2 

OIF •.RI---.-R2 — : 

TEMP « DSQRT (DIF»D1F ♦ BETASq(h>) 

R l » R 1 2 ♦ TEMP 

R 2 * R 1 2 - TEMP 

OIF • D A a 5 ( A 2 - Rl) - DABS(A2 - R 2 1 
IF (OJF .LT. 0.00) GO TO ^20 : . , 

-SHIFT - • -R 2 • — r ! 

GO TO 310 
SHIFT ■ Rl 

go to 310 - - — : : : 1_ 

E 1 G ( l ) * GAMMA ( | ) ♦ SUM 
DO MHO J • 1 »N 

- 1 P 0 S V- (- J 1 _»-J 

IVPOS(J) ■ J 
1 OROI J) a J 

-M. a N - 




i 66* GO TO I/O 

167* **50 DO 4 60 J - I *M 

168* IF <t|G(J> ,GE, E I G ( J* 1 J I- GO. -TO 0*0 - 

I 69* TEMP ■ tJG( J) 

170* E 1 Q ( J) ■ ElG< J+l I 

| 7 I • — — flfcUMl ■ TEHp- ~~ - 

17 2* I TEMP ■ I OHD I J ) 

I 7 J • 1 OKD ( J ) « 1 OKD { J® J ) 

17 4® I OHO I J* t ) ■ I TEMP - 

» 7 B • 4*0 CONTINUE 

176* 470 M » M - 1 

| 7 7 •-- IF (M .NE. Q> 60-T0 v 150 — — — - 

17 8® IF ( N 1 » EQ . 0} GO TO 500 

179® NL » ~ NO I M 1 

180® - J. DO 490 L ■ J s N i 

J8I® NL * NL ♦ NOlM 

182® NV « IORDIL) 

183® NP « I P 0 S V ( N V ) — — — — — _ 

J 8 4 • IF I NP .EQ. L) qO TO 490 

185® LV » I VPOSU.) 

186® I VpOS(NP) » LV — •• • — -- 

187* IPOSV < LV ) » NP 

188® NNP ® INP- 1 I | NO 1 M 

189® r- -DO 400 I 1»N — 

190* j . 1L “ NL ♦ I 

191® INP * NNP ♦ I 

192®- - TEMP o VEC(-IL) — ---. — — 

193* VEC ( IL ) « VEC I INP) 

194® 400 VEC ( INP) « TEMp 

j 95® 490 CONT I NUE — — — — 

196® 500 ESUM - 0.00 

197* ESSQ « 0.00 

198* DO 550 NRR o 1 t N — - — — — - 

19 9® K • N J 

200® 510 K • K - 1 

201 ® IF (K - .EE. 0). GO— TO— 54-0 — 

202® SUM ® 0.00 






203* DO 520 I » K » N 1 

204® :• IE «• (NRR 1 ) * ND I M -♦ — i — ♦- — t— — — — 

205® I J ■ (K-l 1 * ND I M ♦ 1 ♦ 1 

206® 520 SUM ■ SUM ♦ VEc ( 1 L ) • A ( I J ) 

-20-7 •— SUM- n - SUM -♦-- SUM : — : 1 — 1 

208® DO 530 I * K » N j 

209® IL *> (NRR - 1)*N0IM ♦ I ♦ 1 

2 1 0 — IJ ■ • -( K • 1-) *N0 I M I --I — r — — _____ : 

211® 530 VECUL ) = VEC(IL) - SUM*A(lj) 

212® GO TO 510 

2 1 3 * 5 40 ESUM- * - ESUM - f G ( N RR 

214® 550 ESSQ a ESSQ ♦ EIG(NRRJ**2 

215® TEMP o D A B 5 ( 1 2 8 . D 0 * TRACE I 

•DIAGNOSTIC®. THE TEST FOR EQUAL I T y . BETWEEN-nON-INTEGERSMA yNOT-BEmEANINGF-UL*- 

216® IF ( (DA-BSf fRAC-E - ESUM) ♦ TEMP > - TE M P .NE, O.dQ) » N-D ■ 1N0 ♦ 1 

217® TEMP ■ 256.00*EN0RM 

* Q-I-AGNQS-T-4-C-* — T HE- -TE S T— F-OR -E QUAVl T-Y-~bE T-'NEEN---nON^1-N-T54»E-4S---MA y — NO-T--BE-— NE-A-NT-T'LGFTT t-. 

218® IF ( CDABSCEN’ORh * ESSQ) ♦ TEmp) - TEmP #NE. 0*00) |N0 • INO ♦ 2 

219* IF (INO .NE, Oj WRITE (6,86) INO 

220®— 8 6 F 0 R M A T— < ® O®-* — P-R 0 BAB L E — E R R OR — | N— ST-M Q R * • — I-N 0-J-CA-T-O R •-F-.4-5-)- 


22 l * 560 RETURN 7 
22 2 • END 

END OF COMPILATION} 3 DIAGNOSTICS, — - 

.. . = : V - ‘ \.3 

•• ‘ ’’ ' 

it 1 7 6 SEC," ’ ~ " W~ 

92 • • . • ft; 

— 796 7 *i . ~ 

' • * * • • 1 *» • . . • ' . - ■ ' ' 


C OMP i t A T I ON T I ME I 
I /O REQUESTS* 

I/O—WORDS TRANSFEREOi— 

SYMQR . .. ' 

0 1 8 A "09/0 8* I 0 1 93 
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F. MATIN 

Program which solves the set of equations AX = B , and inverts A , 
by means of Gaussian elimination with no pivoting. 

CALL MATIN (A, N, B, NB, DET, PIVRAT) 
where: 

A - N X N matrix of equations (replaced by A 'S 
N - size of A 

B - N X NB matrix of right-hand side of equations (replaced by solu- 
tion matrix X) 

NB - no. of different right-hand sides (if 0, B is ignored by routine) 
DET - returned determinant of A 

PIVRAT - set to ratio of minimum to maximum pivots (-log .^PIVRAT 

is approx, no. of decimal figures lost in the computations). 



J 


Pi^fcl MA TIN 
ri r (ion com 


YCL~ OOP COMPILED 9V. FOR- V S6D1— 0LM-EL9 AUG. JO A T . JJ_L3A.:JL6.«_ 


0101 


SUBROUTINE MATIN t A » N » 8 » N8 » DE T » P I VR A TJ 
L 


mu 


* i ii 1 U U W u = U « 

NX = NX 

4 N 

XX r NX 

4 X 

PIVOT : 

A ( X X_1 

PIVRAT : 

‘oabscpivoti 

IF (AL . 

GT. PIVRAT) AL = PIVRAT 


301 2 M 
25 


OET = OET * P 1 VO T 

IF < P I V R AT .IT. l.D-50) PIVOT 

A_(_K_KJ 

PIVOT = 1 .DO /PIVOT 


r l.D-5 



KI r K - N 
DO 200 I '=' J .NB 

K T r K T 4 N 


B ( X I ) = ( K I ) 4 P I VO T 

DO 1000 J = IfN 




JI : JI ♦ N 
X I : KI 4 N 

I T » - Al l 


IF (NB .EQ„ 0) GO TO 1000 



DO 500 I = I .NB 
J I = J I 4 N 


M 5 ♦ 500 BUI) r BIJI1 * B ( K X I • T 

MR* 1000 CONTINUE 

M 7 ♦ „^.PIVR 

MS* RE TURN 

M 9 * END 


END OF COMPILATION 


NO DIAGNOSTICS, 
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G. STUDNT 

This function returns the value of the student-t statistic for IDF 
degrees of freedom at the 95% level. (Equal-Tails). 

T = STUDNT (IDF) 



IfSTJDNT 
000 COMPILED 


BY FOR-V S7EG 


SEP 71 AT 10:43:59. 


IRAGE USED! CODE { 1 ) 000070J DAT A (0 ) 000052 


1* \ 


FUNCTION 

' z * \ 

C 

RETURNS 

3* \ 


DIMENSI 01 

4* 


OATA T/l 

5* ' “ 


18.2.201*; 

6* 


20.2.074* 

7* : 


IF (IDF 

8 * 


STUDNT r 

9 * 

50 

RETURN 

10* 

100 

IF (IDF 

11 * 


STUDNT r 

12* | 


GO TO 50 

13* ? 

200 

IF (IDF ■ 

14* " 


STUDNT = 

15* 


GO TO 50 

16* 

300 

STUDNT - 

17* 


GO TO 50 

16*. 


END 


( IDF I 


.LT. GO I 
' 1.960 


GO TO 1 00_ 
2 .4 /IDF 

60 TO 200 
2.52/IDF ' 


T ( I OF) 


TNDT' OFf C 0 HP I LATION V 


NO DTAGNCTS' 


COMPILATION TIME: ~ 

1/0 requests: 

I/O WORDS TRANSFEREDS 


0.262 SEC. 

46 
5674 4 
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H. EXAMPLES 

* 

The example is 48 cross sections Q(v) for 1 < E <40 generated 

by means of XSECT (see Sect. II) from a L-J (12,6) potential with 

o 

e = 13.3 cpe, r = 4.92 A (intended to mimic Li-Hg) , y = 9.69925 ppg.» 
m 

g Q = 0.186299, a ± = 0.421559, C g = 0.378xl0 6 cpe-A 6 . (a - 0.25%) 

I. True value of Y- = 0 

d » 

(a) NP = 5 (No Cg) 

(1) Note that the correlation between Y->(C,) and the other 

I o 

parameters is small. 

(2) The "STANDARD ERROR OF ESTIMATE" is 0.262%, which is compatible 
with the precision of the calculated Q(v)'s. 

(b) NP = 6 (C assumed present) 

O 

(1) STANDARD ERROR OF ESTIMATE IS 0.259%, no change from NP = 5, 
indicating Yg is not necessary to the fit. 

(2) The correlation between y. and y, is high. (-0.98) 

1 o 

(3) The determined value of C 0 is not significantly different 

o 

from zero. 

(4) The precision of has decreased markedly. 

2. True value of y^=2 

(a) NP = 5 (C assumed zero) 

O 

(1) STANDARD ERROR OF ESTIMATE = 0.388% (larger than before) 

(2) Precision of C is lower than before. 

D 

(b) NP = 6 .(C D assumed present) 

o 

(1) STANDARD ERROR OF ESTIMATE = 0.260%, lowered from NP = 5. 
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(2) Cg is found to be significantly greater than zero. 

(3) Cg still has comparatively low precision. 


It is to be noted that the presence of Cg , either real or imagined, 
results in an uncertainty in the value of C- . 


TEST RUN-ON 0 * OR I EN 


IAWR- 


DAT ft POINTS 


REDUCED MASS 
INITIAL EPS 
INITIAL RM 

-INITIAL— A 1 

INITIAL GO 
INITIAL CG 


4 969925 + 01 
• 1ROOOO+ 02 
• R700G0+01 
"vR'RO'OO 0+, 0 0~ — 
• 190000+00 
* 378000+06 


WEIGHTING SWITCH 


FIT WITH 5 PARAMETERS CNO C8I 



FIT TO 48 OBSERVATIONS 

THEREAR E 5 P ARAM E-T E R S- IN - THE— MOO EL - — 


§' 

I 


INITIAL PARAMETER VALUES 

— - — -LONER- B OU NO- PARAMETER— VALUE— — — UPPER-BOUND — 

************ **#**************#*****e********************** *•* #****♦.*♦* 


Ml * 


% 


m- 

t*:. 


. 0000000000+000 

•D CO COCOCOG *000 

— rfjosoa oao os +oca 

.0000000000+000 


• 8 3S 8 5 144 57 + 001 
.140 CO COCO 0+0.0 2 
-.'10381-3 094 3 + 003 
.3494053018+004 




.ioooQoaaao+003 

•3CCCC0QCCC+D02 
-5-2000000000 + 005: 
•1QCC0CC0D0+0C5 


a 

fe ■ 


INITIAL SUM OF SQUARES 


.6924430832-001 


ITERATION STOPS ... RELATIVE CHANGE IN EACH PARAMETER IS LESS THAN 


.1CC00DCCCG -CO 2 


1.000000 
— .12 57 35 
.119972 
•0S9570 


l.COOOCC 
-.785783 
-• Si 9235 


1.CCC0C0 

.774383 1 .000000 


“Treee 


-determinant 


.4 S3 932 4 92 2 -02 2 


l-TTG— &F— PI-VO-T-S— 


-^1 G5-tOG-32SS-G £-5~ 







TERM VALUE 95 PERCENT HALFWIDTH 


CG .377715+06 .729954+03 

GO r-2-1-1-93 5+ 00 ; • &27 9 5 8-0 1- 

EPS .139949+02 •102333 + 01 

A1 • 4 35199+00 .3G334G-01 

— - RM- - — ^ 5 r 3529 +0 1 — — — « 687 149+00 


Al+EPS + RM .276225+02 .848760-01 


RM+EPS++2 •888272+03 

: * , .• ’• ' ft . ‘ j 

•351800+02 



4 £ 

FIJ WITH 6 PARAMETERS IWITH 

If; ' . 

C8J 

i ?-, 4 

Vs. 





— ENTRY— TO GASS 5-9 


r IT TO 48 OBSERVATIONS 

-THERE -ARE ~ ~6 — PARAMETERS IN- THE MODEL- 


INITIAL PARAMETER VALUES 


:r.-Bouno 


-V-ALU E UPPER— 3Q-UM C 




.OCCOOCOCOC+OCC .8849899637+0.01 

.0000000000+000 .1399492168+002 

-.000000 CCOC+CC 0 rfrO4-7-T05-4OO+3O 2 

.0000000000+000 .3359161074+004 

.0000000000+000 .GOO 00 00 00 0+0 GO 


.1000000000+003 
.5000000000+002 
-.-20000000 0 C+D 05— 
•1COQ0ODQ0C+005 
•1DOCGOOODO+OC9 


INITIAL SUM OF SQUARES 


.2956066462-003 


ITERATION STOPS ... RELATIVE CHANGE IN EACH PARAMETER IS LESS THAN 


CORRELATION MATRIX 


.loooococca-ocz 


-ROM 1- — — 1-.0&CQ00- — 

ROW 2 .109323 1.000000 

ROW 3 -.067173 -.787775 l.OOOOCC 

-ROW 4 i.11 2 322 — ?5"2 GOT 6 ^-TrttS+S- 

ROW 5 -.096477 -.767335 .964894 

ROW 6 -.980106 -.136736 .092670 


** • v~ - Tt"'' 

.904 T'BB'^'" .-G-tG-O 0-0 tl s* 


.128555 


.113359 i.ocooaa 


OE TERMINAL T 


.129848491 5-027 


RATIO OF PIVOTS 


,1049823447-005 


-SUM- OF- SI.- RESIDUALS .2 327158327 -003- 

MEAN SQ. RESIDUALS .6731329349-005 (FOR 

STANDARD ERROR OF ESTIMATE .2594430555-002 

-SUM-OF-WEI GHTS . 7 - 377 9 5 94-58-00 4- 


42 DEGREES CF FREEDOM) 
CRIT ICAL STUDENT-T VALUE 


o 






NO. PARAMETER VALUE STANDARD ERROR CCNFID. 3 AND HALF-WIDTH 


.1343020206+004 
.3835959553+001 
rl398 6 35876 + 00 2 
.1047772472+003 
.3370098131+004 


.2573184744+001 
.1873918117+000 
• 5 05 37 93 36 6 + C 00 
.1595028352+000 
.6609026839+002 


.5192685754+001 
.3731555717+000 
srlCT. 985 5490 + CCl" 
.3220735173+000 
•1333701S01+GC3 



C6 

60 

EPS 

.375283+06 

.211061+00 

.139864+02 

- .3 6-2751+ B 4- 
.720543-01 
.101986+01 

at 

RM 

.454211+01 

.839337+00 

Al» e~PS*RM 

-276243+02 

.849153-01 

RH*EPS**2 

.888520+03 

.351628+02 

8ETA 

.101129+01 

.148513+01 






•33352 C+OE 


-56-1013+ CL 




' TEST RUN-ON O’BRIEN Q’S -- WI~T H-&8-=~ C 6 — '-t — 

N 0 OF DATA POINTS 48 } 

REDUCED MASS • . 9699Z5+01-- • ~ — 

INITIAL EPS •140000+02 , 

INITIAL RM •470000+01 i 

—I NIT IAL—A1 r44 0000+00— — — 

INITIAL GO .190000+00 , ft 

INITIAL CG • 378000+06 

• : •- - WEIGHTING SWI TCH — .-I——. — . r— • 


FIT WITH 5 PARAMETERS I NO C8I 



r IT TO 48 OBSERVATIONS 

THERE* AR E —5- PARA METER S~ I N- THE— MOEhEt- 


INTTIAL PARAMETER VALUES 

LOUER- SOUND p_araMEI£R— V ALUE URRER--3-OUN-Q 

*****•*•«•*•***• ****************************** ********** ************** 


itiiiuiiiaiiiiiiBiiriyuiii 


.OOGCCOOOOO+OCO 
, CO 00000000*000 
seeeeEeeoeo+cee 
, 0000000000+000 


•03685 14 45 7+0 Cl 
• 1400.000000 + 002 
-. 10 98 13 09 93+-Q- 05- 
•3494053018+004 


iiaiiufl»iHaiii|iiS2aiu:i 


• 1C00000D0D+C02 
.5000000000+002 

-.z&ooeeooee+G es- 

.1000000000+005 


INITIAL SUM OF SQUARES 


.6944308972-001 


— HER ATI G N- -ST G R E-tATI-V E— G KAN-GE-I N — SUM OF — SQ-UA-RE-5 — I-S — LESS — IHAt> 


~^IC ECGC-GG GS-G&; 


CORRELATION MATRIX 


ROW 1 1.CC0000 

ROM 2 —•125349 l.CODOOO 

ROW 3 srl:2C12C ^7855 2:9 2ri 0 D COCO 


•CS9S94 -.619105 

.101897 -.763207 


•774527 1.000000 

• 9 64 7 8C .903471 1.0CC00C 





DETERMINANT .4538018739-G22 RATIO CF PIVOTS .104 8C84199-C05 



SUM OF S^-«— RE-SI DUALS .6 4 81S0S4.S21 - aC3 _ 1 . 

MEAN S9. RESIDUALS .1508152237-004 (FOR 43 DEGREES OF FREEDOM) 

STANDARD ERROR CF ESTIMATE . 3883493578-002 CRITICAL STUDENT-T VALUE 2.017 

SUH-0 F--WEIG HT S .- 733 9488055-0 0 4 — — ^ 


NO. 

PARAMETER VALUE 

STANDARD ERROR CONFID. BAND HALF-WIDTH 

•* ••««** 

1 

2 

. T 

•1346277770+004 

.8842247107+001 

, mm CTTfiii 

•76635638 23 *000 
. 2781089456+000 

.1545437330+001 

.5608357880+000 

ici Q^cnncQ+n^T 


J 

4 

5 

•104 77 7D 37 3+00 3 
.3372315152+004 

.2372077541 + 0 CO 
.9828877413+002 

•4783542557+00C 
.19820 95 972+003 . 


Er 

3ENVALUES OF MOMENT 

MATRIX-FINAL ANALYSIS 



6840167958-002 

5749404385-003 

.2851315924-003 .25 9251968 4 -OCR 

• 1561022 35 4 -C 




t 


TERM VALUE 95 PERCENT H ALFWIDTH 


CG . 377563*06 .106354*04 

GO *212 611 *00 .786564-01 

EPS- - — si 4 016 7*0 2 .1 51635* 01” — — j. 

A1 • 4 35497* 00 •219404-01 ij 

RM .452592*01 .727141*00 [] 

.... ' ‘ ■ ____ -vj 

Al* EPS* RM .276243*02 .126117*00 

R H* EPS ** 2 . 689104*0 3 .522S7& * 02 — — ~4~ 


-V* 


FIT WITH 6 PARAMETERS i WITH C8! 






-IT TO 48 03SERV ATIONS 

-THE-RE—ARE — S RA RAME-TE RS— T N--T HE— MG D EL — 


INITIAL PARAMETER VALUES 

R— '* **. * 

•********«*****+******s****+******+**»***********«*** ***************** 



OOOOCOOOOO+OCO 

QQOOOQQOOQ+QQQ 

D6G-&GO&GOO+-0Oe- 

cooaaooaoo+oca 

occoooooco+coc 


•8842247107+001 
.1401673643+002 
7-70-37 3+ G-G-3- 
•337231 5152+ QC 4 
•OOOOOOOOOO+OCO 




.1000000000+003 
.5000000000+002 
-.GG&EDD&DGC+G 05- 
.1000000000+005 
.1000000000+009 


INITIAL SUM OF SQUARES 


.6485054621-003 


ITERATION STOPS ... RELATIVE CHANGE IN SUM OF S3UARES IS LESS THAN 


«.iaoooocaoc-aa; 



1.000000 
— 1*0843 0 
-.066259 
-.111588 
-- ^ -595 5 96- 
-.980109 


-1-.-000000- 

-.787424 

-.625830 

_ "» c r qc ■? 

© roo jjt 

-.135852 


l.oocoao 

.776739 1.000000 


.091774 


.127835 


-iGeoooe 

.117993 1.000000 


DETERMINA M 


,1264424246-027 


RATIO OF PIVOTS 


•10468183C1-CC3 


SUM OF SG. RESIDUALS .2834185229-003 

-HEAff- SO.- RESIDUALS - — .674-8060070-005 — tF-GR- 

STANOARD ERROR OF ESTIMATE .2597702845-002 

SUM OF HEIGHTS .7339488055-004 


-4-2 DE GRE ES- 0 F - FREE DOM l 

CRITICAL STUDENT-T VALUE 





2 .8828489430+0G1 . 18 780 08742+0 CD . 378 S821599+CCC 

3 .14 010713 20*002 -5057464533*000 .1022514331*001 

4 .104-7-84 2052*003 .-1502 860529*0 00 .-3234572 7 13*000 

5 « 33735208G4*004 -5534937007*002 .1338940363*003 

6 -4127919393*001 .1547427233*001 .3122708120*001 



EIGENVALUES OF MOMENT MATRIX-FINAL ANALYSIS 



.64145 852 2 1-B&2— . 57 44 8 6 0 5&9-6&3 -. 2853 44 11 46-623 ^742674141-664 .-7539-925903 - 925 - 


TERM VALUE 95 PERCENT HALFWIDTH 

**«***«*««*«•*«*.***«******. *********** *********** 


CG 

.375257*06 

GO 

.211848*00 

EPS 

•x'iUrL / *u 2 

A1 

. 435184*00 

RM 

.453093*01 

Al* EPS* RM 

.276261*02 


.354102*04 

.724291-01 

il02261*01 

.490069-01 

.838421+00 


.852789-01 



R H*EP5*~*2 . 8 8942 2*0 3- .-353009*02 


BETA .197094*01 .149093*01 



C8 


739608*06 


566678*06 
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VII. JWKB 


This main program computes JWKB phase shifts ^Tj_ f° r 

reduced energies E = E/^ and values of the orbital quantum number 
such that j£r z. ( JL + ~ ) / ^ CL is a multiple of 0.05 (For details of 
the calculation of > see function PHASHT) . 

The total cross section is given by 


00 


QCE) - 


£ CZJ.+D 

£-0 


(VII-1) 


Letting 


Sg_ - Jr 


(VII-2) 


and 


Since 


i L = ? h 




X 




- 


1. 


(VII-3) 


a - 


k. ^ 

8tt 


, the program computes Q by means of 


a L + aij 


(VII-4) 


where L is a value large enough so that ai= a^-aL can be estimated 
accurately. 

Values of are calculated at every 0.05 in b until one 

of the following is true: 

(1) b > 2 and <0.1 A^/b (This condition 

°2 

forces the remaining cross section AQ to be less than 0.05 A ), 
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3ir < 2 

(2) The Jeff reys-Born long-range phase shift 
and agrees with ^ to within 0.01 radians for 3 successive values 
of 

When either condition is satisfied, this value of JL is taken as L . 

The intermediate values of 'Y! a r e th^n calculated hv interoola— 

l . 

* O 

tion (DNTERP) . The remaining cross section AQ = lil Z\ l, is 


The program could be improved by replacement of DNTERP with a continued- 

is 

, * •* ~ ' “ it 

I Jr 


fraction interpolation program, since the limiting form of 


estimated two different ways: 

oO 


(x) AO, ~ — ^ Jr a2*i \ dlL (QJBVAL) (VII _ 5:) 

L Jl x 

L. . jg 

(The factor of ^ is to correct for the discrepancy 

between and -^iJ0 ) 

(2) V l. is obtained by means of the s-algorithm 


(EPSALG) . 


00 JL-?oo 


Finally, the apparent CA constant is calculated by means of: 

o o 


C 




/ G^. \% 

A ' 8.083' 


(VII-6) 


INPUT : 

Card 1 
Card 2 
Card 3 


Title (20A4) 

NE, MU (15, F10.4) 

(ELIST(I), 1 = 1,NE) (10F8.0) 



VII-3 


where: 

NE - No. of energies (in reduced units) for which cross sections are 
to be calculated 

MU - reduced mass of system in ppg. 


T7T T CT 1 / T ) T _4- V> DnoTOXT 4 t*i imi f o 

tlUXL.' X J. ^ WU CilCl.5) -Li.* JU V_- V— UiiJ. 


NOTE: In addition, a user defined function POTIN/POT is required (q.v.) 


OUTPUT : The values of E , , and a (zero of potential) are in 

reduced units. The quantities EUNIT = e and RUNIT = a are the reducing 
units of energy and distance returned by POTIN (q.v.): they should be 

of the order of e and r^ , respectively. 

* 1 

For J2, corresponding to every 0.05 in b =■ (£ + the 

following are outputted: 


L - JL 
* 

B - b 

YO - y = 1/r 

°° 2 2 

RO - r the outer most zero of 1-b / r - V(r)/E (in reduced units) 
o 


ETATAB(LL) - 


ETAJB - 


1 


3"fi> 




DELQ - & -k-jL re< ^ UCe< ^ un: *-ts) 

IER - Error return from determination of turning point. If = 0, 
normal convergence. If 5 ^ 0 turning point is incorrect (see 


REGFAL) 

. If = 0, normal 
convergence. If 4 0 , did not converge (see PHASHT) . 


ISW - Error return from determination of lx 
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A table of the above quantities is printed out, the units being those 
obtained treating EUNIT and RUNIT as dimensional units. 

Finally, the e-algorithm sequence of diagonal convergents 
(EPSVEC(I) , I = 1,M) are printed out, M being the number of applications 
of the algorithm, and the convergents for J- qq appearing at every 
other element (the first such appearance being the one used). The 
following are then printed out: 



and finally, 

0 'IT 

CROSS - Q - ~ loo 

EPSVEC(J) - 1^ 

C6APP - C* PP 
6 

NOTE: Values of outputted are in reduced units. 


REMARKS : (1) A user-defined potential function- -POTTN-/POT is required (q.v.) 

(2) EUNIT and RUNIT are treated as dimensional units except for 

°2 

actual cross sections which are in A , 

(3) The subroutines TIMSET (in line 54) and TIMGET (in line 93) 
are specific to the UWCC Univac 1108, and set and read an internal clock. 

These functions are not necessary for execution of the program. Line 54 
may be deleted and line 93 replaced with "600 CONTINUE" 
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(4) The convergence tolerances for the calculation, of 71 , 

J8 

convergence to '71 , and determination of the turning point, are 

given by EPS (absolute), TOLJB (relative/absolute), TOLYO (relative), 
in line 8. 

* 

(5) The phase shift calculations generally terminate at b <3 

°2 

and give a cross section which has limits of error +■ 5A . 



JWKB 

T 0 5 / 3 l / 7 .2 w 01 { 53 ? 12 ( , 0 ) JWK B 


1. C PROGRAM which COMPUTES JW«B PHASE SHIFTS AND TOTAL CROSS SECTIONS 

2. C A CONSISTENT SET OF _U NIJS poR T H JS _P R 0 GRAM A R E _ 

3. C ANGSTROMS, PICOPICOGRaMSiAND DEC IPICOSECONOS (GIVING CENTIPICOERGS 

*4 • C _ AS THE UNI T .Q£_.E.N E.R_G 1 1 , _ 

5. REAL MU 

6. DIMENSION E L I S T i 2 0 QJ , E T A J_A BJ..2 0 0 ) ,D I F (2000 U E P_S V_ E C t 2 0) .... . k 

7. COMMON B,RO»E,X 

8 • D A T A E P S / 0 ,005 / , T OL J B / 0.. 0 j / , T O L Y.O / l. £.5.5./ 

9. DATA PI2/1 „57079633/,HX2/l , U2IA2*U/ 

10. . EXTERNAL VE.F.F. 

H. WRITE (6,2} TOLYO ,EPS , TOL JB 

1 2 o . 2 FOR M A T < .1 H i , 5 X , ’ENTRY. TO . J W K B / D N TERP— » / / f 1 5 X , ’ IOL E R A N C E F 0 R Y 0 (RE 

13. ILATIVE) » .EI2.5/5X, ’TOLERANCE FOR ETA ( A BS ) ’ , E I 2 * 5/5 X . * TOLER ANCE FO 

..1**« 2R JB LIMIT? _tE.12.fJM - 

15. C READ IN TITLE OF EXPERIMENT 

... 16 . READ J 5,3 ) J.DJ F.(JJ.., .L*.1.,.2.Q ) ... „ 

17. 3 FORMAT ( 20 A R ) 

1 8 . .... .... .. WR I T E 16, 5.1 (D IF (-U , I.* 1-, 20.) J__ - 

1 9 » 5 FORMAT ( 1 H 1 , l OX , 2QA** > 

-20, . C READ IN NO, OF .ENERGIES NE--AND-REOUCED-MASS_MU - 

21. READ (5,6) NE,MU 

22. 6 .FORm A.T (.1 5 ,.F 1 0 *.H 1 : __ 

23. C READ IN LIST OF ENERGIES 

2*4. READ (5,7) (ELI S.T.U.) .» I a 1 » NE ) 

25. 7 FORMAT ( 1 OF 8 . 0 ) 

26. C POT IN. SHOULD- RETURN THE VALUE OF THE ZERO OF THE POTENT I AL . JL 

27. C C 6 SHOULD BE VALUE OF LONGE*RANGE C6 IN REDUCED UNITS, THE 

.... 28 , C SCALE -OF -WH I.C.H.. ARE....G I V. E N IN E.U.N.I.I.(£NER.G Y ) A.ND—R U.NJJ.JL ENGT.H.) 

29. SIGMA • P0TIN(C6,EUNIT,RUNIT) 

30. .... SIGlNV.a. Is/SIGMA — : 

31 . SIGJN2 « 2,?S161NV 

.32, CONS JB„ ib._.0 ,1 875*P..I Z9-C6. 1 ... . . 

33. ECONS » 2. *MU*EUN I T/HX2 

3 *4 «_ WRITE- (6 ,9.) MU.CA ,SI_GMA , E.UN I-T-, RJUN.I-T- - . 

35. 9 FORMAT ( ////5X , ’PHASE SHIFTS AND CROSS SECTIONS FOR REDUCED MASS’ » 

36. 1 E l-*4. 6 / 5X, ’ LONG-RANGE. C 6 ’ , E lA. 6/5 X , * ZE RO .-OF. -POTeNT I AL ’ „» E 1*4 . 6 / 

37. 2 SX, ’UNITS, OF ENERGY AND LENGTH ’ »2E l *4 *6 ) 

.. 38. DO 1.000- J «. 1, NE i *. 

39 . , E 5= EUSTII) 

..-.SO. WAVfr..®_. SQRT ( E*EC-ON-S l S — .11. 

*M« A ® WAVE«RUN1T 0 

. *42. C JB « CQNSjB*A/E - — 

S3, COEF ** 16«oPI2»RUNlT/WAVE 

*4*4. . WR I TE < 6 , 10 ) - E , WAVE , A . . 

SB. 10 FORMAT (///5X. ’ENERGY* ,F1Q, 5, 5X, ’WAVE NO , » , E 1 5 , 6 1 5 X , ’ A«W A V E *RWN I T » 

_ *4 6 « l o E 1 R « 6 / / 5 X , * FOLLOW I NG TAB LE— XN -REDUC ED -UNI JS EUNl T , RUN IT’/) 

*47, WRITE (6,15) . 

. *4 8, — 1-5 FORMAT. ( 1 *4X . ’L* ,--5x , ’B’ , Q X , ’ Y 0-’-, 1 OX , *R0 A » -9.X , ’ eT A f , 14 X » ’ E T A JB » » - 

*4 9. I 9Xo’DEL9* ,7X,»lER’ ,3X,’ISW’/5X,100(1H«)/I 

50*. 1 JB «* .0 . ; . 

51. C CALCULATE PHASE SHIFTS aT EVERY .05 IN REDUCEO IMPACT PARAMETER 

5 2.. L S T-E P— » A/-20. — — , 

53. IF ILSTEP .LT* 1 ) LSTEP « I 

- _5 <4 « C A L L -T-I-M SE-T-I-0 «-> -- 

55. DO 500 LL a I ,200 




T 





56. NETA a LL 

57 « u —a LSTEP* ( LL-i->— ; — 

58. 8 * ( L + 0 « 5 ) / A 

5 9-. 1 F — ( B . 6 T S-L6MA 4— GO -T 0— i-OQ — — — — x. 

60. XL o SlGINV 

6 1 . XU s .S I g J-N2 — — — 

62. GO TO 200 

63. - -100 XL a 1 » / B - — . 

6 H # XU a SlGINV 

65 * C FIND -TURNING -POINT -Y0—»--l-/R0 — 

66, C IER a 0, NORMAL. IF NOT. 0, ERROR IN ROOT (SEE REGFALI 

6 7 . ZOO C A L fc~R E G F A L < X L , X U-, X , V E F F i TOii Y 0 , g ,20,l£R4 

68. RO o l./X 

69 ■_ — C CALCULATE PHASE- SHI FT— ETA — — 

70. C ISW a 0, NORMAL. IF NOT 0, THEN DID NOT CONVERGE' TO EPS ACCURACY 

^ — 7 1 . — — T 0 L » EPS / A — - — 

72. ETAjB » CJB/B**5 

73. IF (ETAJB -,-L-T .— -1 TOL- « E T A UB *-TOL- — 

7*1. ETATAB(LL) a A*PHA5HT ( TOL , 9 , I SW ) 

75. — • — OELQ « B* ( S I N ( ETAT AB < LL 1 4 4 **2 — - 

76. WRITE (6,20) L , B i X , R 0 , E T A T A B < L L ) » E T A JB , OE L Q , l E R , l S W 

— 77.- 20 FORMAT (10X. 15, 3F»0.6,-3E-I*(. 6.216) — - 

78. IF ( CQEF ♦B*OELQ .GT« 0.1 .OR. B , L T » 2.) GO TO ROQ 

79. WRITE -(6,22-) - 

80. 22 FORMAT ( /5X , • REMA I N I NG CROSS SECTION IS LESS THAN 0.05 SO ANG»» 

-- 8 1 . - —GO- TO 6 00 

82. ROO TOL « TOLJB 

- 83. IF (ETAJB .L-T, l.) -TOL «* ETAJB*TOL - - - 

8M. IF ( ABSiETaTAB(LL) * ETAJB) .LT. TOL ) IJB • I jB ♦ J 

— 85. IF -( I JB , L--T-. — 3— * 0 R -E-T A-J B — .G-T-.-i-.-f — GO— T-0— 5 00 

86. WRITE (6,23) IJB 

87 . 23 FORMAT ( /5X ,* AFTER- ». 18 ,5X ,♦ ASYMPTOTIC-PHASE -SHIFTS,- ETAJB IS LESS 

8 8 o 1 THAN 1 * ) j 

89. GO TO 600 - - 

90. 500 CONTINUE 

•--91 • WRITE (6,25) — 

92. 25 FORMAT (lOX.’MORE THAN 200 PARTIAL WAVES ARE REQUIRED*) 

93. — 600 — T “-TIMGET(Q) i — 

'9 H , ND a 10 ' j 

95. CALL D1FTAb<ETAT-Ab.NETa,DIF-,ND, EPS- 4 — — i— 

96, NM A X = L S T E P * ( NE T A*- |) ♦ I i 

— ■ '9 7. B S T E P * - LST EP/A 1- > — — 

98, 00 a 0.5/A 

-- 99, M ® 0 — i 

100, I EPS a 0 - 4- • 

I 0 I . — — — ■ Q a 0 » : — -4 

102. DO 700 |_L a I ,NMAX I 

-103 ♦ — — — L— °- L L ■» l — j 

104. 0 a ( L + 0 ♦ 5 ) / A 

-105, .. . ' -TOL- » EPS — . 

106, ETAJB » CJR/B*«5 

-107, IF (ETAJB , LT »— t .4— TOL- o ETA JB*TOL — — 

108. ETA = DNTERP<B,B0,BSTEP,ETATA8,NETA»DIF,ND,T0L) 

- 109. DELQ — a B»(SlN(ETA-l-|** 2 - 

110, Q a Q ♦ OELQ 

— 11 l * — IF (LL .LE.-NMAX 20 .0 R-r- IE P S-- .tQ 4-) — GO -TQ-700- 

112, CALL EPSALG ( Q , M , EpSVEC ) 



113* J 9 2 * 2* (M/2) •' H 

„ 1.1 .It - I F ( M ♦ LE . 2 ) Q 0. T 0 7 0 0 — 

11&* IF ( ABS ( EPSVEC ( J 1 • EPSVEe<J*2)) *COEF .GT. 0.11 60 TO 700 

, 1 16 . i eps a j : : : 

117. WRITE (6,21) (EPSVEC(K) ,Kb 1 ,M) 

~.l 1.0 •_ 21 FORM A T — ( / 5 X » * E P- S-LU 0 N _. A L G 0 R !_T H M C ON V E R GES X-I-10F4-U 8)1 

119. 60 TO $00 

120* 70 Q CONT I HUE v 

121 . 800 B * B ♦ l ./A 

- 1 2 2 • _ — Q J B_ s _ Q *. T A ft E T A » Q JBVAL ( C J B , B 1/E T A JB 

123. CROSJB e COEF *Q JB 

— 1 2 8 .» C R O.SS— a_C 0 E ELf JB. : , 

' 125. C 6 APP * HX2*WAVE« (CR0SJB/8.083 ) ft*2*S/MU 

_.l 26* WR ITE ( 6 , 30 IT# CROSS., Q., CROSJB , QJB ,C 6 APP-,ND .. : 

127. 30 FORMAT </>SXi«TlME ( SEC ) « . F 1 0 » 3/ 1 OX , « CROSS SECT I ON AT (.AST PHASE 5 

-1 2 0 » 1HIFT * » F 1 0.2. 10X, »Q» , F 1 0 .3 /l OX, * JB EXTRAP OL A.T £ DC ROSS SECTlO N ( SQ_ A . 

129. 2NG) • .F10.2 , 10X , »QJB 1 .F 1 0. 3 . 1 OX , *C 6 A PP A RE N T • , E 1 S . 6 / 1 OX , * H I GMES T OR 

— 1 30 . _ 30ER -D I FFERENCE . USED J.n I NTERPOLAT.IjON 1 . 1 8 ) 

131. IF (M . EQ « 0) 60 TO 1000 

...132. ; CROSS, s C 0£ F ft EPSVEC (- J ) - - 

133. C6APP n HX2«WAVE* (CROSS/8,083 ) ••2.5/MU 

.13 8. W R 1 1 E _.( 6,80 1- CROSS , EPSVEC (JM C 6 A P P.._ - 

135. 80 FORMAT ( 10X . »EPS1U0N ALGORITHM CROSS SEC T I ON » , F J 0 • 2 , 1 OX , » Q » , F 1 0 . 3 . 

_ J 3 6_» — , l. 1 0 X ..1C 6. A P PA R E N.T » , E 1 8.61... 

137. 1000 continue 

I 38. STOP 

139. END 


ND OF COMPILATION: NO DIAGNOSTICS, 
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A« PHASHT 


This function computes the value of the JWKB phase shift 

1\ 


for 


given E and b = (i + -j) /k , to a given absolute tolerance, by means 
of the formulae 

\ = T, 

■' (V.II-7) 


( >6- > r 0 ) 

Mi. 


where V(r) is the potential, r^ is the largest positive zero of 


F(r) 


F (r) = I - 


Jr 


VO) 


\~ 2 - E 

and y = 1/r (y^ = l/r Q ). The quantities and V are given 

p Yz 


by 


1 Ji /> r * 

Ty = r ^ F Vi - Vl-> 1 ^ 




(VII-8) 


r - 

The integrals are performed by function ADSIMP (q,v.) 

CALL: ETA = PHASHT (EPS, NITER, ISVJ) 

INPUT: EPS - Absolute error tolerance for integrals in Eq. (VII-7,8) 

NITER 

NITER - 2 is the maximum number of points which may be used 


in the integrations. 



VII- 7 


OUTPUT: ISW - = 0, Normal convergence of integrations; if ^ 0, integrals 

did not converge (see ADSIMP) 

PHASHT - value of “3/ (b,E) 

NOTE: Values of B(b), R0(r Q ), and Y0(y Q ) are obtained through the 

COMMON statement of line 10. 


!m phashT 

>•0 s / 3 1/7 2*01 1531.17 .( .0) ...PHASHT ; I v 

| ' l« FUNCTION PHASHT ( EPS ,N ! T£R , I SW ) 

. . 2 . c find? .value of..jw<b_ phase shift_.f.or..^!wp.ac_t„p.arameter.b,. turning.. 

3 • C POINT R0» EPS JS ABS ERROR CRITERION# 2«*NITER IS MAX NO, OF 

R t C - P Q I N TJ . .WH.I C H M AT B E_ U S E D ... .1 N ... J NT E GR A I-I QN (_N LT ER... » l £.*—.?.) 

5» C IF JSW ,NE, 0, FAILED TO CONVERGE TO ACCURACY EPS 

6 . C N 0 T E ... « R E.Q U I RES _ F V N C T I 0 N 5 F . , F 2 , F 3 ., _.W H I CH. F.OR_. V _# .1 /.R... RETURN 

7, C F(Y) m « t "SORT { 1 * ( B*¥ I **2 - V ( R } /£ )) / ( Y * Y ) 

— -0 •— C F 2 t Y.j s SQR.T < 1* ( B« Y ) * * 2* V I R > /EJ-/-I-Y » Y ) . 

V, C p3 I Y ) # (SQRTI l®(B»Y)»*2*»V(RI/E)*SQRTU^lBeY)® 6 2) )/tY«>Y) 

- 1.0 • COMMON B , RO-i e , yq 

11* DATA PJ2/1 ,57079633/ 

■-•.-..-.12 1 E X T E R N A L-- F , F 2 , F 3 _ , 

13# LIMjT « 2*®NJ TER 

1 M i IF. ( B- . L E-» - R 0 I _ 6 0 _„T 0— 1 0 0 — 

15, BINV * 1»/B 

L 16, T 0 L- s -0 *.S * E PS : 

17, E T A 1 * ADS|MP(BINV,YO,F2»TOL* LIMIT • I S W 1 ) 

I B •- E T A 2 * A 0 S I M P ( 0 . , B l N V , F 3 »-T OL » L I M 1 T- »-! 5 W 2 I 

19, 1 5 W * 1SW1 + 1SW2 

.20,.-. . PHASHT « ETA J - * ET.A2 * , 

21. RETURN 

2 2 , UUJ P HAS H .T__ a ... P I 2 » B_*— R.O.* A D S f Mp(0#-,-Y. 0-. F-, EP S...U-M J-T. f 4-5 Wi 

23, RETURN 

iNp OF COMF1 L A T I 0 NI. ... 


NO DIAGNOSTICS, 
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B. VEFF/F/F2/F3 


Four entry-point function which returns values of 

VEFF(y) = 1 - b 2 y 2 - V(r)/E 

F(y) = (1 - Y veff )/y 2 

F2(y) = VEFF /y 2 

F3(y) = ( VEFF ) /y 2 

where y = 1/r , and b and E are obtained via the COMMON statement 
of line 2- 

CALLS Z = VEFF (X) 

Z - F (X) 

Z = F2 (X) 

Z = F3(X) 


ARGUMENT: X = y 




It FUNCTION VEFF(X) 

- 2 t — C 0 MMON B | Rp t E * YO 

3 t U * B * X 

— :—i» t V E F -F- “'-l* — »_-U-«.U— J* — POT-(l t /X)/ E. 

5# RETURN 

- 6 t r — ENT R Y- - F ( X ) — — _ : 

7 t IF (X tUTt I • E - 3 » GO TO 200 

hr- — 8 t U 2 ■ b X * X- ■ - 

9. V s POTU ,/X)/E 

- 10 • WORK - s B*B*U2~ V — 

lit IF (WORK « LT t OtOl) GO TO 100 

— 12* — IF (WORK * 6T * 0 *99999 9 9 * A ND * - WORK — »tT~* — It 0 M— WORK-* 0*9999999 

1 3 1 VEFF ■ (It «■ SORT ( J • *. W0RKII/U2 . 

'•■1 — - 1 <4 ■* : — R ET URN — 

15* 100 VEFF ■ 0 « 5 *(B*a*V/U?)«(i 9 4 0 , 25 • WORK • ( 1 . ♦ 0«5*WQRK>) 

~ 1 6 1 — RETURN T 

17* 200 VEFF *» 0 * 5 *B • B 

■ 18* ■ RETURN- 

■■ 19* ENTRY F 2 ( X 1 

{---2 0 * — — U -• • B * X -- — - — — - 

21* WORK ■ 1, - U*U • POTU./Xl/E 

2 2 * 1 F- ( W 0 R K - u T-* -0 , A N Dr-- W 0 RK— * G-T-r— • 0 , Q l ) - W 0 R K-*-—* 0 R K 

23* VEFF * SQRT I WORK ) / ( X*X ) 

2 <4 * R E T U R N - — — • 

j 25* ENTRY F3(X) 

U--2 6-+ — 1 F ( X *LT *-l . E-3 1--G0 -TO «00 , 

27. U2 b X*X 

28 * BB~w B«B : 

29. OOOA ■ » BB * POT ( l # /X I /(U2*E> 

r 30* — WORK » U2«OOOA- : 

31* IF (WORK *LT* 0.01) GO TO 300 

w— 3 2 * - V E F F * ( S Q R T ( A B S ( 1 . * W 0 R k ) ) * 5 Q RT(A BS<!*» B B • U 2 ) ) -)-/ U 2- — — 

33* RETURN 

‘ 3R *- — 300 WORK a ~DOD A *000 A : 

35* V * BB*BB 

36* VEFF- • -0 « 5«t ( B B*000 A-l — T— 0 * 1 25 * U2 • M V» WORK )- ♦ - 0 * 5*U2 * ( BB* V*> WORK - 

37. 1 OOOA) ) 

j-! — 3 8 * — —RE T U R N 

39. <400 VEFF * 0. 

U — <4 0* RETURN 

• <41* END 


:ND OF COMPILATION! NO DIAGNOSTICS, 

! ’ 

i 




C. POTIN/POT 


This function is a user-defined routine which returns values of 
potential constants (POTIN) and values of V(r) (POT) 

1. POTIN 

CALL: SIGMA = POTIN (C6, EUNIT, RUNIT) 

OUTPUT: C6 - values of C, in reduced units EUNIT = e and RUNIT - a 

6 

EUNIT - e dimensional unit in cpe of energy (of order e) 

o 

RUNIT - a dimensional unit in A of distance (of order r 

m 

POTIN - returns value of 0 (V(a) = 0) in units of RUNIT 

2. POT 

CALL: V = POT(R) 

INPUT: R * r in reduced units 


or a) 


OUTPUT: POT - V(r) in reduced units 
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3. SAMPLE PROGRAM 


The program listing is of a function which reads in a pointwise 

potential (EUNIT = e and RUNIT = r ). The value of V(r) is calculated 

m 

by quadratic interpolation 

V(r) - V -f A r -v By* 2 " 

(VII-9) 

where A(I) and B(I) are calculated by the routine. Extrapolation to 
low r is done by extrapolation of the closest quadratic, extrapolation 

^ 8 , 

r* (vn-io) 

where C8 = C_ is calculated from the last data point and . 

o o 

INPUT : Card 1: NOB (15) 

Card 2: EPS, RM, C6, SIGMA (10F8.3) 

Card 3ff: (X(I), V(I), I = 1,N0B) (10F8.3) 

where 

NOB - no. of points 

EPS - e (cpe) 

RM - r (A) 
m 

C6 - Cg (reduced units) 

SIGMA - 0 (reduced units) 

X(I) - I-th value of r . (reduced units) 

V(I) - I-th value of V (reduced units) 


to large r is done by 


Vtr) ~ 


V* 



REMARKS: (1) C 0 is calculated from C_ = - r^V - r^C, for 

o oD 

V = V (NOB) 

(2) It is not necessary that EPS = e and RM = r 
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r = X(NOB), 


but any scale 


factors may be used. 



POT 

• 05/3 1/72*0 1 15 3 | 2 S | » 0 I POT - : 1 

If FUNCTJ0NP0T1N(C6C0NS,EPS,RM) 

■ r 2 • D 1 MENS.J On — X ( 5Q ) , V < 50 ) , A ( 50 ) , B ( 50 I 

3. READ (5,10) NOB 

— - H t 1 0 — FORmA T ( . J 5 ) 

, 5. READ (5,20) EPS , RM , C 6 , S 1 GMA 

---■*• C600NS * Of — , 

7# 20 FORMAT < ) OF 8 • 3 ) 

-a* RE AO (5,20) (X(l),v(l), !*HNOB) 

9, R 2 « X ( NOB ) *X ( NQB ) 

- — | 0 9 - — RA * R2«R2*R2 

11# C 8 e *R2« ( V { NOB ) ®RA * CA> 

— 1 2 « — R 2 n— l »/X (NOB-1 ) **2 

13. R 6 » R 2 * R 2 * R 2 

18. TEST-* •RA* I R 2* C 8 - C A ) / V I NOB- 1 )--— - 

15, NN b NOB - 1 

1 6 « DO - 90 I - <5 2 , N N 

17, HI m X ( 1 ) » X ( 1 » l ) 

1 8 , — H2 b -X ( 1 + 1 ) «X-1 |4 

19, HH « X ( I * l ) • X ( f •» l ) 

- 20 , D 1 2 s V < l ) -• V ( 1 1 )- — 

21 , D23 8 V ( | + 1 ) « V< I) 

— 22 , HI 2 *-Hl/H2 

23, A ( I ) * ( H l 2*023 ♦ 012/HI 2 )/hH 

— 2*4 , 90 -BID ® ( H 1 2*023 » -D l 2 ) / ( H 1 *hH) — — 

25, NN » ( NOB * 1 1 / 2 

- 2 A , RUP b XI NOB-1 - — 

27, POTfN • SIGMA 

— 2 8 » — R E T U RN ■; , 

29, ENTRY POT(R) 

... 30, — IF (R ,GT, RUP) - 50 TO 500— 

31 , 1 b NN 

32 •- 100 IF ( R ,GT,-X ( l) 4- GO TO 20Q 

33, l » ( !♦! )/2 

—3** , I F- < I — » GT ,-44— 60— TO- 100 

35, I m 2 

3 A • 60" T 0—9 00 

37, 200 J » I ♦ 1 

... 3 8 , j F-M - , E Q # — N-O-B 1—6 0 - T O- 3 0 0 

| 39, IF ( R«X ( 1 ) ) 900,800.200 

[-80 • ——300 1 8 NOB »— 1 

I 8 1, 800 HH » R - X ( !) 

82, POT 1 N b M ( J ) +-HHb(A ( ! )+Hh*B ( ! 44 

83, RETURN 

i- 8 8 * 500 R 2 b 1 • / ( R ♦ R I , 

! 85, RA ■ R2*R2#R2 

-8 A , POT IN s. .R A • ( R2*C6+CA) 

87, RETURN 

.... 8 8 , - END- -r r — 

10 OF— COMP-I (, AT I ON J ... NO DIAGNOSTICS# 



D. QJBVAL 


Function which computes 


tX) 

Sjr 


/ - 2. J8j / ' 

Jr/PW V ciA' 

U 


/ 


2 TR 3tr 

the integral is done by expansion of sin and ^ _ 1 

X 32 

CALL: Q = QJBVAL (CJB.B) 


INPUT: CJB - 3ttAC-/32E 

6 


B - b 


OUTPUT 


: QJBVAL - 


JB 






where 

C^A 

Zjx S 


REMARKS: The value returned is correct to more than 3 significant 

figures for values of B for which -fc ^ £ t 



v \ . 

QJBVAU 


•01 .• 0.1. J.5.3 Ji.7-_U.0J_ QJ.ft.VAU_. 

'■ l « FUNCTION QJBVAUCJBiB) 

C.J 2* JC PE R FOR MS. , I N T EG R A U _F_R.OM._B -TO _ .INF! N m_0 F_ B *.|.S4JL_C-T. A J B »*.#2 

3* C WHERE ETAJB !S J»B PHASE SHjFT FOR UONQwRANGE C6i 

.' Hi... ‘ C ETAJB.. s. CJft/8#«5 i _ C J B *-3 ♦ P ! • C 6 / ( 3 2 *E l : 

'V 5* DIMENSION 5(101 

^.... 6 » DATA S/ | , , • Os 3 3 3 3 3 3 3 3,0 , 0 R RR HR HR R , ®0 « 3.1 7J.6 0.3 i E n 2 • 0 * ! R I 0 9 3 R7.E » 3 

7. I -,R2755S97E'»5,,939&83R5E«,7,*’»IS6A!391E®8,,20H72R06E-J0, 

1. 8* „ 2. ..2 15R9902E" 12/ 

9 , QJBVAU « 0, 

,_1,0 4 T ... s B?a : 

II* Z * CJB/|T#T) 

_ 1 2* Z. « Z*Z 

13* T 0 Z/T 

| r * — 1 D o D a 9 1 — T — : — 

; IS* DO 500 I s 1*10 

L. I 6 • - D E L Q . J?. S.( U_«.D 0 D AAU O f .! e.2 .) 

. 17* QJBVAU » QjBVAL * OEUQ 

_ 1 B • I F ( ABS rOEUQ > -*-L-T.»_l.* E«3*QJB VAU1— G0. T0— 600 

19, 50Q DODA * DODA*T 

- 20* 6.0.0 QJBVAU.. « .Z* QJBVAU 

21. RETURN 

,_22.» END 


!.Q.. OF._C.QMPJ U AT I ON I NO 0JA.GNQ5TJ cS. 
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E. DNTERP 


This function performs finite difference interpolation based on a 
difference table constructed by DIFTAB (q.v.) The appropriate formula 
(Forward, Backward, or Stirling's central difference) is used with maximum 
order such that: 

(1) The order is less than or equal to the size of the difference 
table. 

(2) The size of the next term to be added is less than the prop- 
agated error. 

CALL: VALUE = DNTERP (X, XO, H, N, D, ND, EPS) 

ARGUMENTS: X - abscissa of interpolant 

XO - left-most abscissa of data 
H - stepsize of data 
Y(I) - I-th ordinate of data 
D - difference table from DIFTAB (q.v.) 

ND - highest-order difference available in D (returned by DIFTAB, q.v.) 
EPS - estimate of the absolute error of a Y(I) 

OUTPUT: DNTERP - value of interpolated ordinate' corresponding to X . 

REMARK: The program extrapolates beyond the range of the data by means 

of the nearest forward /backwards difference polynomial. (This could 
be improved by extrapolating the difference table.) 



1 dnterp 

P«05 / 3 i/72 

i i • 

1 i 

i_ 2.« 

! 3* 

0 « 

! 

5 . 

6 8 

7 8 

.... 0 . __ 

9# 

_ l Da 

1 1 • 

- 1 2 8 

i 3 « 

1 M 

IS. 

J 6 . 

17 , 

18a — 

19. 

-20# 

2 1 a 

22. 

23 . 

--2 0. 

25, 

26. 

27, 

-28 0 

29. 

30. 

31 . 

32. 

33 . 

-30. 

3 5 . 

36, 

37, 

38® 

39 , 

00. 

01 a 

02. 

03 . 

00 . 

0 5. 

— 0 6 « 

07. 

08. 

09. 

50 , 

5 1 ® 

5 2 a 

53. 


-i-N—yj-S- 


0-U-5J-? 18 ( .0) 0 N T ERP- 

FUNCTION DN T E RP I X e XO*H»Y »N»OsND»EPS) 

C- — - — P £ R F-0 R MS DIFFER E N C E— IN T E R P-0 LA -T4 ON -UP -TO-ORO £R— NO 
C D - DIFFERENCE TABLE FROM DlFTAB 

C E P S— - ES T-I-M A-T-E- 0 F_ A B S 0 LUT&- ERROR 

DIMeNSION Y (J ) , D I l ) 

M = 0 , — — 

DODA ® 1. 

TOL EPS-- — . 

THETA = (X«XO)/H 

I-e -THETA -♦ --0-.-5 

IF (I .LT. 0) I <= 0 

J F - <1- , GE . N-) — l—O-N—w — J 


I 


THETA o THETA 

1— a -I — ♦ -J — 

dnterp ■ » y ( n 

— JF < I « GT . — N~._ NO) — GO— TQ -2000- — 

IF (I « L T » ND) GO TO 1000 

— USE ST IRL I NG * S CENTRAL D I FFERENCE - FORMUL A- 

DO 000 J » 1 ,N0 

- T 0 L -• 2 . * T 0 L — — 

L ' « J / 2 

— K - »- -M- ♦ I -•?- L 

IF (ABSID(K) 1 .LE. TOL) GO TO 500 

— DODA - a DODA / J - - — 

IF ( 2 * L »NE» J) GO TO 300 

— TERM ® THETA*DODA -- 

DODA 0 DODA* ( THETA-L ) 

— GOTO -350 

DOOA = OODA* ( THETA + L > ■ : 

— DNTeRP * D N T E R P— ♦ - 0 • 5« doD A * ( D I K >—♦ — d I K * l I-)—-- 

GO TO 000 ) 

— DNTERP-** DNTERP-* -TERM *D <K> !- 

M a M * N ■» J i 

— R E T-U R N L _L 

USE NEWTON’S FORWARD DIFFERENCE FORMULA i 

-00 1 500 J at -l-.NO - 

TOL * 2.eT0L ! 

— K * -H * I- — r — — i- 

IF (AQS(D«K)» .LE. TOL) GO TO 1600 

— DOOA «- D 0 D A * I T H E-T-A *• J *4 )/J i- 

DNTeRP • DNTERP + DODA*0<K) j 

- M -= M + -N « J — — — j- 

return ‘ 

— - U 5 E N E W T ON * S B A C K W A RD S — - D I F F E RE NC E— FOR HUL A— i- 

00 2500 J a 1,ND ; 

— TOL--“ 2 . » T OL — L 

K « M * 1 „ J 

-IF lABS(DlK)-) .LE. TOL) GO TO— 2600 - 

DODA « DODA* I THETA + U« l ) /J 

— D N T E R P— *— ON T E R P— ♦— D 0 0 A • D I K > — — 

M b M ♦ N - U 

— R ETU R N 

END 


300 


350 
000 
.500 . 
C 

1000 


1500 

1600 

.c . .... 

2000 


2500 

-2600 


END OF COMP J LAT I ON I 


NO diagnostics. 
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F. DIFTAB 


This program constructs a finite difference table from N equally- 
spaced data points, with ordinates (1 1 i £ N), Successively higher 
order differences are computed until no element of the current column 
in the difference table is larger than the propagated error, 

CALL DIFTAB (Y, N, D, ND, EPS) 

INPUT: Y (I) - y I-th ordinate 

N - no. data points 

EPS - estimate of absolute errors in Y(I) 

OUTPUT: D - difference table stored in packed-mode: i. e . , the first- 

order differences are in D(l) ... D(N-l), the second-order 
differences in D(N) ... d| 2N-2^ etc. 

ND - set to highest-order difference calculated. 


: 


diftab 

• 05/ 3 1/7 z rn Q 1 J J, 3 1 20 .1,0 ) 
1 

c 

c 

c 


JJIFTAB 


2 

3 

.... H 

5 

6 

7 

8 
9 

.JO 
1 1 
I 2 
13 
_ 1 H 
15 
U 
17 
. 1 8 

19 

20 
21 
22 


SUBROUTINE OIFTAB<y,N,0, NO, EPS) 

C A L CUt ATE S.OI FEE REN C E. _T A B L E OF 1 ^ _ 

NOTE • 0 MUST BE DIMENSIONED AT LEAST NO*(N«U 
ND-«- U P.0 N .. E NT RY , M AX .0 R D E R..D ! F F RE N C E -TO- BE - C A tC UVAT E 0,-UP 0 N 

c return is Actual' highest order 

C E PS * ABSO t u T £ ERROR ESTIMATE. 0 F_ ¥ A S-T0-S TO R.._D IF F.ERENCES 

DIMENSION YU ) ,01 1 ) 

M A X 0 -.' N D - _ — . 

ERRLVL * 2 # • E P S 

T o FOROIF ( Y,M,0) 

ND : . 

K m 1 

.100 .... I F (NO .6 E,- M AXO ,0R, -M . » L E .1 , 0 R. _T ,tE, ERR LVlJ. G 0. TO 200 

ND m NO ♦ J 

ER Rl y L.B.2«.tERRtVL_ : _ _ 

J * K ♦ M 

- _ — T _m._F.0R01.JM 01 K.) , M , D-U ) ) 

K • J 

GO TO 1 00 ; ; 

200 RETURN 

E ND 


NO OF. compilation: 


. __N 0. — D-l~A-G.N0 STicS. 


I 





G. FORD IF 


Function which calculates the vector of forward-differences 
corresponding to a vector of ordinates. 

CALL: DMAX = FORDIF (Y, N, D) 

INPUT: Y(I) - I-th ordinate 

N - No. of data 

OUTPUT: D(I) - I-th forward difference, Y(I+1) - Y(I) 

FORDIF - max | D (I) | 1 < I < N-l 


! FORD I F 

)aOB/i 1/7 2 « 01-| 5312 1- _ ( , 0 -F-0 R P l F- — — 

i l* FUNCTION fordif (Y»N,DI 

} - - 2*. C C A ic u u TES FORWARD --D f F F E R E N C E S--D ( I )— s — -Y4-I-* 1-1 — •— Y-'-t I 4 

| * 3t c returns VALUE of maximum difference 

} H* DIMENS I ON Y C 1 -| »D< 1 1 

5 » N s N <* | 

7* DO 100 I s hN 

8. - D ( I > *■ y(|*|)-.-4|.M 

9. 100 IF (T 9 L 7 o ABStomi ) T ■ ABSIDMM 

; 4 0 • FORD I F — P-T - 

i lit RETURN 

• — t 2 • — —END * 

|N D 0 F- COMF1 L A T I 0 N 1 -r-- N 0 01 AGNOSTICS. — 
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HL EPSALG 

Subroutine which extrapolates a sequence of partial sums to a limit 

A 

by means of the e-algorithm , Only the diagonal of the matrix of con- 
* 

See P. Wynn, "Five Lectures on the Numerical Application of Continued 
Fractions", Orientation Lecture Series No. 5, Mathematics Research Center, 
University of Wisconsin, Madison, Wisconsin, 53706. 



CALL EPSALG (S, M, X) 

INPUT: S - current value of partial sum 

M - Number of times EPSALG has been called previously for this 
series 

OUTPUT: M - will be replaced by M + 1 

X(l) ro. X(M) - will contain the e-diagonal of convergents ,= 
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REMARKS: The approximations to the limit of the series will be in every 

other X(I) . The best value for the limit will be found in X(l) 
if M is odd or X(2) if M is even. 



EPSALS 
•05/3 1 /72« 

0115 3 
C 

122. j ,01 EPS ALG ... . 


1. 

2 , 

subroutine; epsalgisiM.xj 

extrapolates sum of series s by epsilon algorithm ... 

3 t 

C 

C 

s is value of current PARTIAL SUM. at exit 
T HE. E p S I L ON. . 0 J.A GQNAL. XJ..1J.. . ...X (.M l 

M IS LENGTH OF 

St 
6 • 

c 

BEST ESTIMATE Of SUM IS IN X ( J ) » WHERE J m 
DIMENSION XI 1.) . . . .. _ . 

2* M0D(M,2I 

7 * 
8, 


K f M 
M o M ♦ 1 


9, 

10, ... 


X < M | m 0. 

A 1 a S 


i 1 1* 

100 

IF <K ,LE» 0) GO TO 200 


12, _ . 


AO b X ( K 1 | ♦ i ./ < A i »X ( If » > 


13. 

. 1*». .... 


X ( K ♦ 1 1 *» A 1 
A 1 m A 0 


1 5 , 


K « K « 1 


! 1 6 , 


GO TO 100 


17. 

_ 18* . 

200 

X (1) m A 1 
RETURN ... ... . 


19. 


END 



ND OF COMPILATIONS NO DIAGNOSTICS, 
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I. REGFAL 

This subroutine finds a root of the equation F(JC ) = 0 , given 
an interval which contains the root, by a strategy using regula falsi 
and the method of bisection. Barring multiple roots or excessive 
rounding-error, the routine guarantees convergence. The algorithm pro- 
ceeds as follows: 

(1) Bisect the original interval until a sign change is found. 

(If none is found, error return.) 

(2) Apply regula falsi (secant method) until convergence is 
attained. 

(a) If the change in interval length at any step is less than 
10% of the previous interval length, apply the method of bisection once. 

(b) If convergence is not attained, or the root is lost, 
error return. 

CALL REGFAL (XL, XU, X, F, EPS, LIM1, LIM2, IER) 

INPUT: XL, XU - Endpoints of an interval containing the root 

F - name of function which stores the value of F(OC ) in y in 
the call Y = F(X) 

EPS - relative error criterion for root. Termination occurs when 
relative change in root is less than EPS or the root itself 
is less than EPS 

LIM1 - Limit on number of applications of bisection to find a sign 

change; if LIM1 = 0, Step (1) above is skipped 

\ 

LIM2 - limit on number of applications of regula falsi in step (2) 
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OUTPUT: X - estimate of root 

IER - 0 Normal 

1 No root or a multiple root in interval (to tolerance EPS) 

2 No sign change found with LIM1 iterations 

3 No convergence to EPS tolerance with LIM2 iterations 

4 Root was lost in regula falsi iterations (no longer a sign 
chang e ) 

REMARK: The convergence testing could be improved slightly by adding 

line 12.5 HO = EPS*H 
and replacing line 50 with 

IF (ABS(DELX) . LT . HO) GO TO 1000 



regfal 

"06/3 1/-72 
1 • 

2 . 

3 , 

M 9 

5. 

_ 6 • — 

7. 

... e. : 

9. 

— 10 » 

1 1 • 

12. 

13. 

... | M . 

» IS. 

_ 16 , 

17. 

. 18 , 

19, 

..20, - 
21 , 

— 2 2 * 

23, 

2 M » 

25, 

26, 

2 7. 

— 2 8 , 

29, 

.. 30, — — 

31, 

- . 3 2, 

33, 

— 3 M . 

35, 

36, 

• 37, 

38. 

39, 

0 • 

Ml, 

M2. — 

M3. 

MM. — - 

M 5 . 

— M 6 * 

. M 7 , 

- - M 8 . 

M 9 , 

-60 . 

51 , 

-52. 

63, 

— & H , 

65 , 


, 044 - 53-1 


2 M (.01 RfGFAl — 

SUBROUTINE REGFAL ( XL , X U . X , F , EP S » L I H 1 . L I M 2 » 1 ER 1 

F I NOS-ROOT-OF -THE-EQTN F LX-I- - - 0—1 N -I NTER V AL-! XL , XU I — — 

8 Y METHOD OF REGULA FALSI ( GUARANTEED CONVERGENCE ) 

EPS .» RELAT 1-V E - £ R R 0 R— T OLE RAN C E— FOR . -ROOT 

L I M I - ALLOWABLE ITER. JO F I ND SIGN CHANGE IN AN INTERVAL 

L I M 2 —*• ALLOWABLE -I T-ERAT JONS -TO FIND-ROOT , — 

1ER a 0, NORMAL. »l NO OR MULT ROOT TO TOLERANCE EPS, ■ 2 NO 
SIGN- CHANGE « l TH- L I M 1_ -J-TE R A T-l 0 N S , a-3 N OCON VER G EN C E»l T H L1M* -|Y. 
® M, LOST ROOT IN REG. F ALS I 1 TER AT I PNS 

FA. a F( XL) — -- - 

X » X L 

H a -XU .-XL'- : 

IF (LIMI ,EQ, 0 ) GO TO 300 

J *» -. 1— — -1 — 

HH e EPS *H 

DO 2 00 I 1-tL-I'Ml 

IF (H c LE » HH) GO JO 250 

x *_ X L 

DO 1 00 K » 1 * J 

FB F ( X ♦N) - i — — 

IF (FA*FB , LE » 0*1 GO TO 350 

X a - X -H : t 

H » 0 ♦ 5 «H 

j .« 2 * J 

I ER « 2 

RETURN — 

IER a 1 

return 

FB a F(x+H) 

B » X * H 

F A - m F ( A ) — — r — 

HH o ABS(H) 

XP -b -A — ■ — ■ — t 

DO 800 1 a |,UH2 ! 

X- «- A • FA • H /-( F B * Fa I i~:_ 

J B 1 

FX -o - F ( XT — i — * 

IF ( F X A 1 500,1000,600 i 

FB-a FX : r-T- 

B a X 

GO TO 700 — * i — 

FA a FX - - ! ' 

IF ( FA*FB .GT, 0») GO TO 900 

H - a - B «* -A 

IF ( J ,EQ. 2) GO TO 800 

OELX » X- ■* - XP 

XP m X 

IF ( A 8 S ( DE L X ) — .UT-,-EP S • A B S(-X)—. 0 R •— A B S< X)—«LT«E P S-I— G 0— T 0-1000- 

I F (ABS(H) ,LT, 0 * 9 * HH ) GO TO 800 i 

X -a — 0 « 5* ( A* B4 : * t~ 

J a 2 j 

GO— T-0--MOO 

HH P ABS { H ) I 




800 




54 • I ER » 3 

.1. .6.7 j_ RETURN 

58# 900 I ER s H 

59# RETURN ■ _ 

60# 1000 1 ER 8 Q ; 

1 6 1 t R E T u R N 

' 62 t END 



NO OF COMPILATION! j NO OIAQnOSTICS. 
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,J. ADSIMP 


This function returns the value of the integral 

.xu 


FIX) dX. 




obtained by means of adaptive application of Simpson's rule. The program 
is a slightly modified version of Alg. 182 , Comm. ACM 315 (1963). 

CALL: VALUE = ADSIMP (XL, XU, F, ERRTOL , MAX, IER) 

INPUT: XL, XU - Endpoints of the interval of integration 

F - name of function which stores the value of the integrand in Y 
from the call Y = F(X) 

ERRTOL - Error bound on absolute error of the integral 
MAX - maximum number of integration points which may be used 
OUTPUT: IER - if 0, normal convergence. If ^ 0, then more th^n MA2( 

points are needed for convergence. 

ADSIMP - Returned value of integral 
REMARK: All of the dimensions in lines 6-8 could be safely reduced to 

15. (which necessitates a 15 in line 41) 



ADS J MP 

>05/3 L/J7-2 * 0 1 I 5 34 2 a f » 0 I — AOS I MP . ... 

I* function ads IMP ( XL , XU ,F .errtol »max , I er ) 

i 2. C — INTEGRATES F-( XI-FROM XL-TO-XU BT_.ADAPT.IVE SIMPSONS RULE 

j 3. C ERRTOL » ALLOWABLE ABSOLUTE ERROR TOLERANCE 

• R » __C— MAX -r .LIMIT— ON N 0 _F-U N C-T I 0 N—C ALLS r 

; So C • IF I E R • NE » 0 MORE THAN MAX POINTS ARE REQUIREO 

; 6 B __ ; D I M E N S I 0 N- D x ( 3 0 ) ».E P S E c 3 0 ) . X 2 I 3 0 I , X 3 L3 Q.) , F 2 ( 3 0) ,p3 (30) , F*M 30 ). 

i 7. IFMPI30) »FBP(30) , E S T 2 ( 3 0 ) * E S T 3 ( 3 0 ) »PVAL(3Q»3) 

>-• Bo — INTEGER -RTRNI 30) , : — 

9. EPS a ERRTOU 

1 0 o A _® X L__ : : 

1 I o B 0 XU 

.. ■ 1 2 o LVL ■- 0 .. 

I 3 o EST m l oO 

— 1R. DA -o-.B m ..A ; .. 

IS, F A a F ( A ) 

— 1 6 * — FM D R . Q * F l-C A -kB-l-*.0. 54 — ; — 

I 7 o FB a F < B ) 

... | 8, KOUNT-a ■ 3 — ; : 

19. 100 L VL » L VL +1 

•—■2 0 « D X ( L VL ) * Q« 3 3 3 3 3-3 3 3 3 * D A — — - — 

2 1 o SX » DX ( LVL I *Oo 1 66666667 I 

— 22. FI o -R « 0*F ( A - ♦ -OoS*DX (LVL )-) — 

23. X2 ( L VL ) ■ A ♦ DX(LVL) 

— -2Rt— - - F 2 (LVL) p F(X2lLVL-l) 

25. X 3 I L V L ) o X2ILVU * DX(LVL) 

— 26# F 3 ( L VL ) a F ( X 3 ( L VL-) ) ~ 

27. EPSp(LVL) » EPS 

— 28. FR ( L VL ) -®— R , Q»F-( X^3 (VVL-) — t— Q , 5*0 X4L-VL- )-) — 1 

29. FMP(LVL) « FM 

— 30. EST l « SX* ( FA ♦ F I -* F 2 ( L VL ) ) *— 

31 . FBPlLVL ) ■ FB 

- 3 2 . E S T 2 ( L V L ) « S X * ( F 2 ( L V L -I- * F 3 < L V L > — *— F M-) 1 , 

33. E 5 Y 3 ( L V L I '■ SX*(F3(LVL) * FR ( L VL ) * FB) j 

—3 R * SUM .a EST4 - + EST2-( L-VV)-*-EST3 (-LVL ) j—- ■ — : 

3 S o KOUNT a KOUNT ♦ R j 

- 36 « I F ( ABS < EST«SUM )- • 5T •- EPSP-( U-VL ) - ) GO — T0--500-J 

379 ROO LVL « LVL w I 

. -38* : 1 a. RTRN (LVL)- — — 

39c P V A L (LVL o I) = SUM 

... .R0« GOTO (600, 7 0 0 ,800 ) ,-J — — 

Rio 500 IF (LVL . GE 9 30 ,0R* KOUNT qGgo MAX) GO TO ROO 

- -R 2 o RTRN (LVL) m — I 

R 3 9 FM a F l i 

• — R R • - FB a F2 ( L Vb ) U-, 

H 5 o EST a ESTI 

— R 6 * 550 DA o OX (LVL I : i 

R 7 o EPS a Oo577«EPSPUVL) > 

— - R 0 9 GO TO 100 — * — 

R 9 . 600 RTRN(LVL) » 2 

50 o FA a F2ILVLT 1 — -r 

51 . FM a FMP(LVL) 

—52. F 8— □— F 3 ( L VLr) - — — 

5 3 t EST a EST2(LVL) 

-SR * A - P — X 2 4 L VW4 

55, 60 TO 550 ‘ 



66 » 700 RTRN (L VLI 8 3 

_.§7 i FA..* F3 ( LVuJ__._ 

1 5a. FM * FatLVU 

.. 59.» F 8 . b_ FBP I LVl.) 

60. EST * EST3 (LVU) 

_1.6 I • ..... A * ...X 3 ( L V L > I 

. 62. GO TO 550 

_63 * 800 SUH PVAU ( LV.L »X)_*. P.V A L-U. V L . 2 ) P.V A UL V.L...3 ).. r 

6*». JF (LVL .GT. II 60 TO 800 

.65. ADS IMP. a SUM : 

66. |ER * 0 

6 7 . I F...I K 0 U N T . . G T_»_MAX_) !. E R-P—l 

68. RETURN 

.6 9 »._ END 

ID OF.. COMP. J L AT. ION | ...NO DIAGNOSTICS. 


i 





VII-21 


K. SAMPLE OUTPUT 


The problem shown is for the Na-Hg system with the pointwise potential 
taken from U„ Buck and H. Pauly, J. Chem. Phys. 54, 1929 (1971). 

o 

The units of energy and length were e = 8.79 cpe and r = 4.72 A , 

m 

& 

resp. The output shown is for E = 1.65. 

jJ* 

NOTES: (1) The last partial wave calculated was at £ *» 455, (b =3), 

°2 

when SQ < 0.1 A /b, 

(2) The convergents of the e-algorithm after 5 extrapolations were: 


- 175.4996 (best), 175.4999, 175.0902. 

(3) The calculation took 6 seconds on an Univac 1108 (1.5 ys add time), 
Since every seventh "77 was calculated, without interpolation the 
execution would have taken 40 seconds. 

°2 

(4) The e-algorithm extrapolated cross section was 696.6 A , with 
Cg PP - 0.6690 x 10 5 cpe - A 6 . 





PHASE SHI FT S AND CROSS SECTIONS FOR REDUCE D MA SS. »3H2513»Q2 
LONG»RANG£ CA •SQOOqO+QQ 

ZERO OF POTENT I A L .fi lOBRQ + QQ ' 

UNITS OF ENERGY AND LENGTH »879000*01 ,«*7200Q«<01 




ENERGY 1 e 6500 0 

wave no. 

( 298886*02 

AsWAV£*RUNlT ,191079+03 




FOLLOWING TABLE IN 

REOUCED UNITS 

eunit.runit 





L B 

YQ 

RO 

ETA ETAU8 OELQ 

IER 

t * a • <# a * • « 

isw 

ios*9««a** 



0 

9003599 

1 .922920 

. 702780 

*•,856023*02 

,225136*19 

.175033902 

0 

0 

7 

9 053 1 63 

1*921865 

a 703302 

*•6 799 108*02 

,296975*08 

s 1 16608901 

0 

0 


. 102783 

1 .91 8988 

.709728 

- a 697999*02 

.109763+07 

,90881 1-01 

0 

0 

21 

c 152902 

1.919359 

.707037 

-.551 189*02 

,153195*06 

,199912*00 

Q 

0 

28 

e 20202 1 

1.908079 

0 7 1 0 S 90 

-.960106*02 

,379171+05 

,162609*00 

0 

0 

.35 

9 25 J690 

J.90Q376 

.7 19099 

-.379170*02 

, 129782 + 05 

,195012*01 

0 

0 

<♦2 

.301260 

1.391112 

.713899 

9 ,293999*02 

.507399*09. 

,233086*00 

0 

0 

*♦9 

.350879 

1.380208 

.729528 

-,217882*02 

,236738+09 

,192985-01 

0 

0 

56 

.900998 

L* 3 6 8 1 9 2_ 

.730892 

. • ,_l 9.7.2 60*02 

. 122195*09 

,276979*00 

0 

0 

63 

.9501 18 

1.355286 

.737852 

-.81 6297*01 

,681937+03 

.908993*00 

0 

0 

70 

1 9997 37 1.391338 

.795529 

_5_u2.0 9 5 0 8*01 

.903969*03 

.379505*00 

0 

0 

77 

.599356 

1 .326383 

,753930 

,397989*01 

,251693+03 

,605077*01 

0 

0 

89 

.598975 

1.310977 

.763081 

,855836*01 

.163309*03 

.397806*00 

0 

0 

91 

.698595 

1,293687 

,772989 

.131965*02 

,109695+03 

, 199877*00 

0 

0 

9_8 

• 698 2 1.9.. 

1 ,276065 

,7.8 3 659 

*. 17 2.909*02: 

. ,758775*02 

•6972J 1 *00 

0 

0 

105 

.797833 

1 .257697 

.795136 

.208903*02 

,538309*02 

,623590*00 

0 

0 

112 

.797952.. 

23833 L 

.807538 

23 9 39 3*02 

.390920+02 

,689169*00 

0 

0 

U’ 

.897072 

1.218999 

• 820682 

.265923*02 

,288705+02 

,82525J*00 

0 

0 

126 

.896691 

1 .198939 

.839922 

__j 286279*02 

,217191+02 

,107989*00 

0 

0 

133 

.996310 

1 . 176721 

.8998 1 9 

,301672*02 

,165916*02 

.851979*00 

0 

0 

190 

.995929 

L* 15 9.330 

.866303 

,311727*02 

l 28503 + 02 

,577692-01 

0 

0 

197 

1 .095599 

1.130917 

.889236 

,316020*02 

,100771+02 

.357779*01 

0 

0 

159 

1.095168 

1 . 105898 

,909289 

U.113 2L*0 2 

.799193+01 

.286879-03 

0 

0 

161 

I . 199787 

l .079500 

,926355 

,306323*02 

,690371*01 

*570906*00 

0 

0 

| 68 

1.199906 

1.050019 

.952368 

.291 102*02 

,517958*01 

.657322*00 

0 

0 

175 

1.299026 

1.017129 

,983165 

,268976*02 

,922580+01 

, 121839*01. 

0 

0 

182 

L.293695 

.978690 1.021779 

,2 3.7 23 3*02 

.397520*01 

,126025*01 

0 

0 

189 

1.393269 

.929633 

1 .Q7S693 

.197253*02 

,287909*01: 

,792292*00 

0 

0 

196 

L« 3 9 2 8 8 3 

,869297 

1.157077 

fLi.9.9 9 5 9*02 

.290199*01 

.599686*00 

0 

0 

203 

1.992503 

.790788 

1.269561 

,103836*02 

,201591*01 

,966929*00 

0 

0 

210 

1.992 1.22 

.,7.25005- 

1 .379301 

•-7.Q7 9 7 2*01 

.170229+01 

.755222*00 

0 

0 

217 

1.591791 

.689871 

1.960130 

,527661*01 

,199593*01 

.110085*01 

0 

0 

229 

L. 5.9 1.3.60 

•.65 9 22 8 . 

L. 5 2 8 5 19 

,905952*0 l 

.123371*01 

.1009 1 7*01 

0 

0 

231 

1.690980 

•628806 

1 ,590316 

.313950*01 

.105813+01 

,826092-09 

0 

0 

238 1.690599 

.609719 1.653675 

,237289*01. 

.91 1701+00 

,817165*00 

0 

0 

295 

1.790218 

.589066 

1,712136 

, 185593*01 

,788927*00 

,160253*01 

0 

0 

252 

L._7 8 9 837 

.565683 

1,767773 

l19.779Q*QI 

,685968+QO 

.177999*01 

0 

0 

259 

1.839957 

.598997 

1.823165 

.117991*01 

,597871+00 

.156592*01 

0 

0 


266 

1.889076 

.533001 

1 

.676168 

,991833*00 

,523369*00 

. 132352*01 

0 

273 

1.938695 

.518827 

1 

.927929 

. 855278*00 

.959735*00 

. 110938*01 

0 

280 

1.988319 

.505919 

I 

.978579 

.739179*00 

,905163*00 

,892986*00 

0 

287 

2.037939 

.992703 

2 

•02962] 

,632512*00 

.358183*00 

.7 12228*00 

0 

299 

2.087553 

.980683 

2 

• 080371 

,598260*00 

.317591*00 

.567089*00 

0 

301 

2. 137172 

.969277 

2 

.130939 

.977267*00 

,282395*00 

.950956*00 

0 

308 

2,186791 

.958930 

2 

.181358 

,917195*00 

,251778*00 

.358956*00 

0 

315 

2.236911 

.998098 

2 

•23)656 

.365977*00 

.225059*00 

,286906*00: 

0 

322 

2.286030 

.938290 

2 

.281859 

,322221*00 

,201672*00 

.229299*00 

0 

329 

2.335699 

.928822 

2 

.331968 

,289637*00 

.181191*00 

.189175*00 

0 

336 

2,385268 

.919813 

2 

,362019 

,252239*00 

,163068*00: 

• 198570*00 

0 

393 

2.939888 

.911189 

2 

• 932001 

.229201*00 

.197116*00 

. 120355*00 

0 

350 

2.989507 

.902911 

2 

.981939 

,199832*00 

,133001*00 

,978999*01 

0 

357 

2.539126 

.399970 

2 

,531835 

. 178609*00 

,120980*00 

.799813*01 

0 

369 

2.583795 

.387392 

2 

.581696 

. 160093*00 

,109397*00 

.656165*01 

0 

371 

2,633365 _ 

.380008 

2 

.631527 

.193761*00 

.999260-01 

.590507-01 

0 

378 

2.682989 

.372999 

2 

•681332 

. 129991*00 

*905659-01 

.997026-01 

0 

~3B5 

..2.7 32603 

.36615J 

2 

, 731.115 

:1 6 80 l_* 00 

.826366-01 

,.371101-01 


392 

2,782222 

.359599 

2 

.780879 

.105632*00 

,755259-01 

.3D9289-D1 

0 

399 

2.831892 

.353279 

2 

.830626 

,957 176- 0 1 

.691369*01 

.258658*01 

0 

906 

2,881961 

.397179 

2 

.880358 

,869193-01 

,633857-01 

.217121*01 

0 

913 

2.931080 

.391288 

2 

.930078 

,790669-01 

.581991-01 

.182855*01 

0 

920 

2.980699 

.335599 

2 

.979788 

,720631-01 

,535136-01 

.159523*01 

0 

927 

3.030319 

.330089 

3 

,029987 

,657892-01 

,992735-01 

.130970*01 

0 

939 

3.079938 

.329762 

3 

.079J79 

,601676-01 

,959303-01 

. 1 1 1369*01 

0 

99 1 

3. 129557 

.319605 

3 

. 128863 

.551199-01 

.919912*01 

.999672-02 

0 

998 

3.179176 

.319610 

3 

. 178590 

.505678-01 

,367688-01. 

.812256-02 

0 

955 

3.228796 

.309769 

3 

.228217 

.96961 2-01 

,358800-01 

,696981*02 

0 

gaining 

CROSS SECTION 

IS LESS 

Than 0.05 

SQ ANG 


1 



S I LON ALGORIT HM C ON V E R : G E S , ; 

175.9996 -59.9791 175.9999 102.8919 175.0902 


i g ( SEC) 5.937 ; ; 

CROSS SECTION AT LAST PHASE SHIFT 699.92 0 175.090 

JB EXTRAPOLATED CROSS SECTION JL5.Q_.ANG ) 696,51 QJB 175. 99 0 C 6 A PPARENT *66 8939 

HIGHEST order difference used in interpolation 10 

EPSILON^ ALGORITHM CROSS SECT I ON __ 696*55 5 17_5,.500 


C 6 APP AR ENT 


.66903L*05 
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